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Adaptive  HF  Propagation  Path  Utilization 


l.  INTRODUCTION 

IIP  propagation  paths  arc  time-spread  and  frequency-spread  channels,  and 
are  characterized  by  severe  variability  in  the  time  domain  of  all  their  properties, 
inclusive  of  path  losses. 

Improvements  over  a  present-day  IIP  link's  performance  in  terms  of  circuit 
reliability,  data  rate,  and  error  rate  can  be  achieved  only  through  the  use  of  adap¬ 
tive  schemes  that  in  principle  appear  to  be  able  to  cope  with  the  channel  variabil¬ 
ity,  if  the  necessary  penalty  in  terms  of  equipment  complexity  is  willingly 
accepted. 

Real-time  oblique  ionospheric  sounding  between  the  two  terminals  of  the  link 
and  the  channel  probing  simultaneously  performed  between  them  are  the  data- 
gathering  operations  that  provide  the  inputs  on  which  to  base  the  adaptive  control 
of  the  link's  performance  parameters. 

For  the  sake  of  illustration,  we  have  assumed  in  this  report  that  the  link  has 
1125  to  3375  carriers,  respectively  for  a  midlatitude  and  a  transauroral  path, 
available  in  the  HP  band  between  3  MHz  and  30  MIIz. 

A  sounding  scan  lasts  100  to  160  seconds  and  is  repeated  every  300  to  480 
seconds.  The  master  station  of  the  link,  where  the  sounding  transmitter  is 
located,  also  generates  the  waveform  for  channel  probing  and  includes  a  complete 

(Received  for  publication 
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terminal  for  two-way  digital  communications.  During  the  pauses  of  the  emissions, 
measurements  of  noise  and  interference  levels  are  performed  at  both  the  master 
and  the  slave  station  of  the  link,  for  use  by  decision-making  microprocessors  and 
control  units.  At  each  terminal,  the  transmitting  and  the  receiving  facility  could 
have  separate  units  for  sounding/probing  and  for  communicating,  or  these  func¬ 
tions  could  be  performed  by  the  same  equipment  in  different  modes  of  operation. 

In  the  latter  case,  the  equipment  at  the  two  terminals  could  be  identical  and  the 
assignment  of  the  master  and  slave  roles  would  be  dictated  by  operational  require¬ 
ments. 

By  processing  the  data  obtained  by  sounding  and  probing,  it  will  be  possible 
to  select  automatically  the  group  of  frequencies  to  be  used  for  communicating.  At 
each  sounding  cycle,  Information  about  the  frequency  selection  and  about  the  wave¬ 
form  to  be  employed  is  exchanged  between  terminals  and  used  locally  to  achieve 
adaptivity.  During  the  next  sounding  scan  (they  are  performed  at  a  rate  of  one 
every  5  to  8  minutes)  the  gr?up  of  frequencies  that  were  selected  for  communi¬ 
cating  are  excluded  from  the  sounding  frequency  plan.  Instead,  Information  on  the 
changing  status  of  the  group  of  communicating  frequencies  Is  obtained  from  meas¬ 
urements  performed  on  the  coded  waveform  that  is  part  of  the  communications  bit 
stream. 


2.  TECHNICAL  DISCUSSION 

2.1  General 

During  the  course  of  this  study  we  have  evaluated  two  basic  configurations 
for  the  link: 

1.  a  two-way,  high-data-rate,  high-capacity,  adaptive  exchange  between  the 
two  terminals  of  the  HF  link,  with  adaptivity  based  on  path  sounding  and  channel 
probing,  with  feedback  acknowledgement; 

2.  a  one-way  link  for  relatively  reliable  information  transfer,  although 
characterized  by  lack  of  feedback  acknowledgement  and  lack  of  adaptivity. 

In  this  report.  Configuration  1  is  Illustrated  in  more  detail.  Configuration  2  is 
treated  in  Section  2.  7,  as  an  application  of  M-ary  codes  with  a  large  value  of  M. 

2.2  Pertinent  Properties  of  HF  Ionospheric  Propagation  Paths 

2.  2.  1  GENERAL  REMARKS 

The  HF  propagation  paths  of  interest  in  this  study  are  single-hop  paths  estab¬ 
lished  in  midlatitude  or  across  the  auroral  zone.  The  properties  of  interest  are: 
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1.  path  losses; 

2.  noise  and  interference  spectral  density; 

3.  multipath  spread; 

4.  Doppler  spread. 

Of  these  properties,  it  would  be  relevant  to  know  accurately  the  time  variability 
and  the  related  statistics.  Unfortunately,  this  is  known  only  in  particular  cases, 
so  that  a  reliable  experimental  investigation  on  the  properties  above  is  thus  far  an 
unfulfilled  requirement.  In  general,  we  can  say  that  these  ionospheric  channels 
exhibit  time  fadings  that  are  important  in  determining  the  design  of  the  signal  and, 
in  addition,  show  long-term  variations  due  to  large-scale  fluctuations  of  the  med¬ 
ium.  Such  slow  effects  have  a  time  constant  significantly  greater  than  5  to  10 
minutes,  an  interval  of  time  selected  (as  will  be  seer*  in  the  following  sections)  as 
the  basic  sounding/probing  periodicity.  Adaptive  approaches  to  the  communica¬ 
tion  problem  arc  required  to  circumvent  this  long-term  variability  in  propagation 
conditions. 

In  this  study,  we  will  make  the  usual  distinction  between  path  sounding  and 
channel  probing,  with  the  former  devoted  to  the  measurement  of  path  losses  and 
of  noise  and  interference  levels,  and  with  the  latter  devoted  to  measurement  of 
such  parameters  as  mult:,  th  spread  and  Doppler  spread.  The  following  criteria 
were  adopted  in  our  study: 

1.  The  link  is  assumed  reciprocal,  except  for  the  noise  and  interference 
levels  at  each  terminal.  Therefore,  the  decision  on  the  frequencies  to  be  used 
(this  decision  is  the  output  of  the  sounding  operation)  is  based  on  the  measurement 
of  Interference  and  noise  both  at  the  master  and  at  the  slave  station,  and  on  the 
one-way  measurement  of  the  path  losses  between  the  two. 

2.  Processing  of  multipath  spread  and  Doppler  spread  (the  basic  operation 
of  the  channel  probing  activity)  is  performed  al  the  slave  station  and  the  results 
are  transmitted  back  to  the  master  station,  for  use  in  the  final  selection  of  the 
frequencies  to  be  used  in  communicating. 

3.  Channel  probing  is  to  be  undertaken  only  at  the  best  frequencies  put  in 
evidence  by  the  path  sounding,  in  order  to  shorten  the  overall  cycle  sounding/ 
probing. 

2.  2.  2  HF  IONOSPHERIC  PATHS  CHARACTERIZATION  AND  MEASUREMENT 

2.  2. 2.  i  General  Considerations 

Due  to  the  propagation  mechanism  of  radio  waves  in  the  ionosphere,  various 
mathematical  models  can  be  formulated  to  identify  its  measureable  parameters. 
These  models  (or  characterization)  are  also  useful  in  selecting  measurement 
techniques  for  the  channel  parameters.  The  select'on  of  channel  parameters  itself 
depends  upon  the  particular  application  and  upon  the  receiver  structure.  In 
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Section  2,  2.2.  G,  we  have  summarized  tne  general  principles  involved  in  randomly 
modeling  time-variant  channels.  These  models  are  the  multipath  transmission 
models  and  the  composite  filter  models.  The  two  types  of  models  may  be  distin¬ 
guished  initially.  A  "multipath  transmission"  view  is  suggested  by  identifiable 
mechanisms  or  "paths"  of  propagation  between  input  and  output  terminals  The 
"composite-filter"  view  is  a  view  whereby  channel  output  is  related  to  input  by 
means  of  a  linear  mathematical  operation  involving  a  suitably  defined  system 
function. 

The  HF  ionospheric  channel,  by  virtue  of  energy  propagating  over  many  pa.hs 
of  different  propagation  delays,  exhibits  time  dispersion  (or  frequency  selective 
fading),  and  by  virtue  of  time  variations  of  individual  phase-path  lengths,  exhibits 
frequency  dispersion  (or  time-selective  fading). 

For  most  communication  problems,  it  is  usually  valid  to  assume  that  the 
channel  is  linear,  with  nonlinear  effects,  if  any,  being  relatively  negligible.  This 
suggests  modeling  the  channel  as  a  time-variant  linear  (TVI.)  filter  in  terms  of  a 
time-variant  transfer  function. 

2.  2.  2.  2  The  Basic  Signal  Models  Employed  in  the  I1F  Channel 
Characterization  (Multipath  View) 

These  can  be  described  as  linear  combinations  of  the  responses  of  three  types 
of  propagation  paths;  namely, 

1.  the  ground-wave  propagation  as  modeled  by  a  steady  specular  path. 

2.  a  random  scatter  mechanism  of  propagation  ns  generally  modeled  by  a 
diffuse  path. 

3.  long-term  variabilities  of  ionospheric  paths,  often  as  simulated  by  a 
quasi -specular  model  which  is  characterized  by  a  log-normal  envelope  and 
Gaussian  phase  probability  density  functions.  In  quasi -specular  type  of  propaga¬ 
tion,  the  observed  fluctuations  in  envelope  and  phase  are  caused  by  fluctuating 
attenuation  and  electric  path  length  over  the  traversed  path. 

The  output  of  a  given  medium  may  actually  consist  of  a  linear  combination  of 
Ron-diffuse  and  diffuse  components,  and  different  combinations  may  be  necessary 
to  describe  the  output  of  the  same  medium  at  different  times.  For  example,  the 
sum  of  a  few  quasi-specular  paths,  each  of  which  fluctuates  slowly  relative  to  the 
fluctuations  of  the  resultant  signal,  is  modeled  by  a  sample  function  of  a  Gaussian 
process  plus  a  specular  component.  Such  a  combination  may  be  encountered  over 
HF  links. 

The  representation  of  a  channel  output  in  terms  of  a  linear  combination  of 
non-diffuse  and  diffuse  signal  models  with  stationary  statistics  tiansfornis  the 
characterization  of  the  channel  into  the  specification  of; 
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1.  A  multipath  structure, 

2.  The  relative  intensities,  and  average  delay  and  Doppler  differences  of 
the  various  distinguishable  paths, 

3.  The  statistical  characteristics  of  individual  path  delays  and  delay  spread, 
and  path  Doppler  shifts  and  Doppler  spread. 

4.  The  constitution  and  properties  of  the  fine  structure  of  each  distinguishable 
path  (that  is,  its  possible  decomposition  into  the  sum  of  non-diffuse  (specular  or 
quasi -specular)  and  diffuse  components)  including  the  relative  characteristics  of 
the  components,  and  the  statistical  properties  of  their  parameters. 

When  modeled  by  a  "composite  filter",  the  variable  multipath  nature  of  the 
propagation  process  over  HF  channels  causes  a  number  of  effects  of  great  impor¬ 
tance  in  the  evaluation  of  tire  signal  transmission  performance  of  the  "equivalent" 
filter. 

First,  the  group  delay  differences  among  the  various  paths  cause  the  overall 
channel -composite  filter  attenuation  and  delay  characteristics  to  vary  with  fre¬ 
quency.  If  the  delay  spread  of  the  significant  paths  is  not  a  very  small  fraction  of 
the  reciprocal  of  bandwidth  occupied  by  the  signal,  the  various  components  of 
the  signal  will  experience  non-uniform  attenuation  and  delay,  which  results  in  sig¬ 
nal  waveform  distortion.  The  characteristics  of  this  distortion  will  be  random 
for  a  randomly  time-variant  linear  channel. 

Second,  the  fluctuations  in  the  relative  characteristics  of  the  various  paths 
cause  each  frequency  component  in  the  signal  to  acquire  a  modulation  of  envelope 
and  phase  (and,  lienee,  frequency).  This  multipath-induced  modulation  causes  the 
signal  level  to  fluctuate  up  and  down,  which  raises  the  possibility  of  signal  outages 
or  dropouts,  each  resulting  from  a  "fade"  or  drop  of  the  received-signal  strength 
below  the  threshold  of  acceptable  performance  in  the  presence  of  independent 
additive  disturbances.  In  addition,  a  non-zero  delay  spread  among  the  various 
paths  limits  the  bandwidth  over  which  the  fluctuations  experienced  by  signal  com¬ 
ponents  at  different  frequencies  will  maintain  the  necessary  degree  of  mutual 
coherence  to  keep  the  resultant  signal  distortion  below  tolerable  bounds. 

Finally,  the  fluctuation  rate  induced  by  the  channel  in  a  particular  parameter 
of  the  carrier  sets  a  non-zero  limit  on  the  frequency  content  in  the  baseband 
spectrum  that  can  be  used  to  modulate  that  particular  carrier  parameter  if  inter¬ 
ference  between  the  channel  fluctuations  and  the  desired  signal  baseband  waveform 
is  to  remain  negligible. 

If  the  "composite"  filter  model  is  used,  the  time,  frequency  and/or  statistical 
behavior  of  a  channel-characteristic  system  function  offers  the  basis  for  the 
definition  of  the  gross  parameters  for  describing  the  gross  aspects  of  the  channel 
response.  The  gross  transmission  parameters  of  principal  importance  m  charac¬ 
terizing  the  effects  of  a  randomly  time-variant  medium  upon  broad  classes  of 
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signals  are  discussed  In  Section  2. 2. 2. 7.  In  the  next  section,  we  will  first  discuss 
the  system  functions  for  fading  dispersive  linear  channels  (or  randomly  time- 
variant  channels),  then  we  will  outline  in  Section  2.4,  the  measurement  techniques 
of  the  gross  transmission  parameters  summarised  in  Section  2, 2. 2. 7. 

2, 2. 2. 3  System  Functions  for  Fading  Dispersive  Linear  Chamtele 

The  system  (Unctions  for  radio  channels  have  been  extensively  studied  in  the 
1  2  S  4  &  fi 

past  '  *  '  '  '  ”,  Below  we  will  briefly  summarise  the  results  In  a  way  as  to 
Identify  important  channel  parameters. 

Radio  channels  can  generally  be  represented  by  randomly  time-variant  linear 
(RTVL)  filters.  System  functions  can  then  be  defined  for  representing  RTVL  chan¬ 
nels  in  the  r.ame  manner  as  Is  done  with  other  types  of  linear  filters.  There  are 
four  ways  n  which  the  output  of  th'e  filter  y(t)  or  its  Fourier  transforms  Y(ju)  can 
be  expressed  in  terms  of  the  input  x(t)  or  X(jw)i  namely 
y<t)  in  terms  of  an  operation  on  x(t) 
yft)  In  terms  of  an  operation  or.  X(ju) 

YQu)  In  terms  of  an  operation  on  x(t) 

Y(Jw)  in  terms  of  an  operation  on  X(Ju). 

The  required  expressions  are  given  in  Table  1,  where  h(t,  {)  la  also  called  a 
delay-spread  function  and  H($u,  t)  Is  the  Fourier  transform  of  h(t,  I),  in  the  treat¬ 
ment  of  s  channel  as  a  signal-distorting  filter  with  randomly  time -variant  charac¬ 
teristic  parameters,  extremely  useful  analytical  and  structural  models  can  be 
deduced  (as  for  s  time -variant  linear  filter)  by  examining  the  consequences  of  the 
dependence  of  h(t,  f )  upon  S  and/or  of  H(]w,  t)  upon  u,  holding  t  fixed.  Therefore, 

In  the  subsequent  discussion,  t  will  be  dropped.  Based  upon  these  two  approaches 
(time  and  frequency  domain  approximation)  two  different  representations  of  the 
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Table  1.  System  Parameters 


inPutx<t>  _  y(t) 

- •»  TVL  Channel  - •»  output 

X(1u)  -  Y(ju) 

l.  Time-smear,  or  multipath-spread  function,  h(t,  f ) 


y(t)  *  /  x(t -f)  h(t,  ?)  d5 


=  t  responses  to  variously  delayed  and  weighted  replicao  of  x(t). 
2,  Frequency-smear,  or  Doppler-spread  function  M(u,  C) 


▼00 

Y(jw)  *  f  X(jw  -  j€>  M(w,  €)  d<€/2») 


-  £  responses  to  variously  frequency-shifted  and  weighted 
replicas  of  X(jw). 

3.  Time-dependent,  frequency  selective  function,  H(jw,  t) 


y(t)  »  /  X(ju)  iHju,  t)  e3*4’1  d(w/2») 


=  E  responses  to  X(ju)  e3*^  time  components. 

4.  Frequency -dependent,  time-selective  function  m(t,  u>) 


+00 

Y(ju)  J  x(t)  m<t,  u)  e'3wt  dt 


E  responses  to  x(t)  e’3ut  frequency  components 


v,  vv.„ .. 


2.  2.  2. 4  Sampling  Model 

Restriction*  of  the  frequency  range  of  Interest  in  the  analysis  of  a  filter  re¬ 
sponse  can  be  interpreted  as  equivalent  to  replacing  the  actual  filter  by  the  same 
filter  in  cascade  with  a  rectangular  filter,  as  shown  in  Figure  1.  The  frequency- 
response  function  H(ju)  of  the  filter  is  thus  replaced  by  H(ju)prect  ^  (u),  where 

prect  fl  ^  represents  the  rectangular  amplitude  characteristic  shown  in  Figure  1, 

The  resulting  amplitude  characteristic  13  of  course  strictly  physically  unrealizable 
and  hence  the  corresponding  models  may  not  be  strictly  physically  realizable. 
However,  in  practice,  this  need  not  cause  any  difficulty  because  of  approximations 
in  the  necessary  number  of  paths  in  the  model. 
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Figure  l.  Truncation  of  a  Filter  Transfer  Function  by  Restriction 
of  the  input  Spectrum  to  the  Range  Occupied  by  the  Input  Signal 


We  now  observe  that  since  the  model  frequency -response  function,  H(ju) 
Prt>ct  p  ^  is  /tM'°  f°r  M  *  !l„/2,  its  inverse  Fourier  transform  (by  the  sam- 

plmg  theorem!  can  be  expressed  as 


sm 


2  if 


^mud.'l^*  X  ** model  (  11  Q 


f(M( 

_  ri2ir  ) 

\  2  ! 

i'  )\ 

n  -<*? 


Hi 


(1) 


This  suggests  the  parallel-path  representation  shown  in  Figure  2.  For  a 
finite  number  of  paths  all  with  positive  delays,  the  filter  can  also  be  modeled  by  a 
tapped  delay  line  preceded  by  a  rectangular  filter  as  shown  in  Figure  3.  The  path 
gains,  and  hence  the  tap  gains,  are  determined  from 


model 


(t)  = 


2x 

2  » 


h(T) 


(ty 


(2) 


by  taking  sample  values  spaced  uniformly  2»/flx  seconds  apart. 
2.  2. 2.  5  Paired-Echo  Model 


In  practical  applications,  one  is  usually  concerned  with  how  well  the  actual 
filter  characteristics  can  be  approximated  by  a  given  number  of  paths,  or  "taps", 
in  the  parallel-path  or  tapped  delay-line  model.  To  minimise  the  number  of  taps 
required,  Baghdady 1  proposed  the  paired-echo  model  which  approximates  the 
system  function  in  frequency  domain.  This  approximation  determines  the  number 
of  necessary  paths  and  the  explicit  interrelations  among  the  tap  gains.  This  is  in 
contrast  to  the  sampling  approach  where  the  approximation  is  made  in  the  time 
domain  with  no  explicit  quantitative  indication  of  how  well  the  result  approximates 
the  frequency-domain  characteristics.  The  paired-echo  model,  discussed  below, 
is  of  practical  value  in  computations  of  distortion  only  when  the  number  of  parallel 
signal  paths  in  the  resulting  model  is  not  too  large.  The  method  also  becomes 
impractical  if  the  nonlinearity  of  the  phase  characteristic  is  quite  severe. 

In  the  paired-echo  approach,  the  starting  point  is  H(jw)  which  is  the  low -pass 
analog  of  the  bandpass  filter  relative  to  the  reference  frequency  of  the  signal. 

For  convenience,  lei 


H(ju)  =  A(w)  e'Ww) 


(3) 


where  A(w)  and  Mu)  are  both  ceal  functions  of  w,  known  respectively  as  the  ampli¬ 
tude  and  phase  characteristics  of  the  filter.  As  is  well-known,  the  filter  does  not 
introduce  signal  distortion  if  A(u)  is  constant  and  Mu)  is  linear  over  the  frequency 
range  occupied  by  the  input-signal  spectrum.  Distortion  results  from  departures 
in  A(u)  and/or  Mu)  from  these  conditions.  These  deviations  A(w)  -  A(o)  and 
Mu)  ~  u>Tg(j<0)  may  be  expressed  as  the  sums  of  approximately  chosen  and  weighted 
trigonometric  terms  within  the  frequency  range  occupied  by  the  signal  spectrum  in 
Fourier  series.  Thus,  if  W  denotes  the  width  of  the  band  occupied  by  the  input 
signal,  then  one  can  write 
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Figure  2.  Parallel-Path  Representation  of  h. 


A(u)  -  A(o)  -  Y,  an  e 


jn(2tr/W)u 


*,u)-ured(o.=  l  el”'2’/W>“  ■  £  Cm.in 

m  =  -oe  m  =  l  L  J 


where  we  have  set  |bm|  A  Cm  and  Mbm  +  */2.  The  Fourier  coefficients 

a  and  b  are  obtained  in  the  usual  manner  with  W  as  the  basic  averaging  interval, 
n  m 

Using  the  Fourier  expansions  of  Eqs.  (4)  and  (5),  we  can  write  from  Eq.  (3) 

r  1  *  -jCsin(m(2*/W)u>i3  1 

H(jw)  -  A(o)  +  T  a  IT' 

^  n  m*l 


-jwTgd<l>) 


for  |u|  £  j 


Each  term  in  the  product  within  braces  can  in  turn  be  expanded  in  an  infinite 


series  as 


-iC  sinlm(2if/W>w *-**  I  “  , 

•  m  ’  I  WJ* 


where  J.  (C  )  is  the  Bessel  function  of  the  first  kind,  order  K,  and  argument  C  . 
k  ru  •** 

From  Eqs,  (G)  and  (7)  and  collecting  terms  having  identical  exponentials,  we  obtain 


Mut  (o)  *  j[K(2ff/W)u.^oKl 

H(3W)  -  e  ^  y  |«KI  e 


where  a,,  is  a  complex  coefficient.  Various  mathematical  models  from  Eqs.  (G) 

K 

and  (8)  can  be  obtained.  In  Figure  4,  we  have  shown  a  model,  dei  ivcd  from 
Eq,  (8),  which  connects  the  input  and  output  terminals  by  a  multiplicity  of  parallel 
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Figure  4.  Block  Diagram  Representation  of  H(ju) 


branches,  each  providing  a  separate  transmission  path  for  the  signal  between 

input  and  output.  The  individual  paths  pass  the  signal  without  distortion;  only  pure 

delays  and  constant  attenuation  and  carrier-reference  phase  shifts  are  introduced. 

It  is  also  Important  to  note  the  symmetry  of  paths  about  K  -  0;  that  is,  the  paths 

for  +K  and  -K  have  conjugate  complex  transmission  “gains",  but  the  delay  of  the 

+K  path  exceeds  r^fo)  by  *he  same  amount  that  the  delay  of  the  -K  path  is  less 

than  r  j(o).  Symmetrical  path  pairs  about  the  K  -  o  are  thus  said  to  contribute 
8°  17 

“paired  echoes'  *  of  the  input  signal  m  the  output.  Each  trigonometi  ic  term  in  the 
expansion  for  A(u)  contributes,  in  the  absence  of  any  distortion  from  4(u>),  only 
one  pair  of  echoes,  a  precursory  echo  (corresponding  to  (+)ve  n  in  Eq.  (41)  and  a 
successory  echo  (for  (-)ve  n  in  Eq.  (4)).  But  in  the  expansion  of  4(w)  in  Eq.  (6), 
each  term  gives  rise  to  an  infinite  number  of  uniformly  spaced  precursory  and 


7.  Bello,  P.A.  and  Esposito,  R.  (1970)  Measurement  techniques  for  time-varying 
dispersive  channels,  Alta  Frequenza,  No.  II,  Vol.  XXXIX,  pp.  980-996. 
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successory  echoes.  From  Eq.  (8)  and  Figure  4,  the  quantity  t^(o)  *  K  (2»/W) 
will  become  negative  for  K  >  W  Tgd(o)/2x,  unless  rgd(o)  si°Pe  of  $(w)  at  u>  =  o) 
is  infinite.  The  representation  In  Figure  4  therefore  Is  a  physically  unrealizable 
structure.  Therefore,  Fourier  expansions  of  A(u)  -  A(o)  and  $(u)  -  u  Tgd(o)  that 
converge  rapidly  are  of  course  most  desirable  from  the  viewpoint  of  truncation  of 
the  series  to  fewer  terms.  Filters  with  very  smooth  and  low  ripple  in  A(u>)  and 
^(u)  may  also  be  represented  by  a  much  smaller  number  of  pure-delay,  all-pass 
paths.  When  the  number  of  parallel  paths  in  Figure  4,  or  in  the  model  for  Figure  5, 
is  finite  and  the  corresponding  delays  are  all  positive,  the  signal-distorting  filter 
model  can  also  be  redrawn  in  terms  of  a  tapped  delay  line  as  shown  in  Figure  6. 
Such  a  model  may  be  physically  conceived  to  consist  of  a  lossless  delay  line  with 
taps  whose  number  and  spacing  are  chosen  equa.'  to  the  number  and  delay  differ¬ 
ences  of  the  parallel-path  model. 


2. 2. 2. 6  Statistical  Characterization  of  Random  Time-Variant 
Linear  (KTVL)  Transmission  Channels 

The  randomness  of  the  channel  results  from  the  variation  with  time  of  the 
properties  of  its  active  propagation  mechanisms  and  therefore  is  reflected  in  the 
characterization  of  time -dependent  system  functions  h(t,  ?)  and  t).  The  sta¬ 
tistical  characteristics  of  the  response  of  a  RTVL  channel  depends  upon  the  time 
duration  in  which  the  observations  are  made  of  the  channel  response  for  the  sta¬ 
tistical  characterization  of  the  channel -induced  "luctuations.  The  selection  of 
averaging  time  intervals,  dependent  on  fading  rate,  is  normally  aimed  at  estab¬ 
lishing  time  intervals  over  which  the  generally  nonstationary  statistical  process  of 
interest  can  be  modeled  approximately  by  a  statistically  stationary  process.  If 
the  nominal  fluctuation  rate  or  fading  bandwidth  is  denoted  by  Bfad,  then  the 
"instantaneous"  averaging  time,  due  to  detection  circuits,  is  taken  as  less  than 
one  tenth  of  the  reciprocal  of  Bfad.  The  time  interval,  which  is  long  enough  to 
include  a  "typical"  pattern  of  the  instantaneous  channel-induced  fluctuations,  is 
normally  greater  than  1000  times  the  reciprocal  of  Bfad.  Traditionally,  this  time 
interval  is  chosen  on  the  basis  of  the  time  interval  mat  yields  received  envelope 
statistics  in  response  to  a  test  sinusoidal  carrier  with  a  best  fit  to  the  Rayleigh 
distribution  function.  This  time  interval  is  important  in  the  sense  that  H(ju,  t)  can 
be  treated  as  wide-sense  stationary. 

A  combined  time-shift  and  frequency-shift  correlation  function  ol  H(ju,  t)  is 
defined  as 


Rjjjj(S7,  r)  =  tune  and  frequency  average  of  {H(ju,  t)  H  (-}w  -  JH,  t  ♦  t)}  ,  < y ) 


A  double  Fourier  transformation  of  R^H  (fl,  t)  over  fl  and  t  yields 

S(Td-fd>  =  /  /  RHH(n*T)  e3nTd  e'j2,rfdT  d  (S?)  dT  (10) 

-os  -» 

If  the  variables  and  f^  are  Interpreted  as  time-delay  and  Doppler-shift  variables, 
then  SiT^,  f^)  is  called  the  channel-scattering  function.  Since  f^)  and  RHH  (fl,  t) 
Rrh  (0,  T)  ®re  Fourier  transforms  of  each  other  with 

n  —  Td  and  r  — *  fd  , 

the  spread  of  R^  (fl,  r)  in  the  fl  variable  is  inversely  proportional  to  the  spread 
of  SlTj,  fj)  in  the  delay  variable  r^.  The  spread  of  R^H  (fl,  r)  in  the  r  variable  is 
Inversely  proportional  to  the  spread  of  S(r^,  f^)  in  the  Doppler-shift  variable,  f^. 

The  scattering  function  in  Eq.  (10)  can  also  be  obtained  from  the  tap-gain 
correlation  function  which,  for  WSSUS  (wide-sense  stationary  uncorrelated  scatters) 
channel,  is  defined  as 

R(C.  r)  =  <h*(t.  €)  h  (t  +  r,  «)>  .  (11) 

where  h(t,  f)  dC  denotes  the  complex  gain  associated  with  scatters  in  the  interval 
(€,  €  +  d?).  It  is  called  the  input  delay  spread  function.  The  Fourier  transform  of 
the  time-varying  impulse  response  h(t,  5)  with  respect  to  5  is  H(jju,  t)  which  is  also 
called  the  time-variant  transfer  function. 

From  the  Fourier  transform  of  the  R(€,  t)  with  respect  to  t,  one  obtains 
S(Td,  fd)  and  if  one  Fourier  transforms  R(€,  t)  with  respect  to  the  multipath  vari  ¬ 
able  f,  one  obtains  R(fl,  t),  defined  as 


R(U,r)  -  J  R(f,  r)  exp  <-j2*fH)  d?  .  (12) 

R(fl,  t)  is  called  the  time-frequency  correlation  function.  This  is  the  cross- 
correlation  function  between  received  earners  spaced  c?  Uz  apart  at  transmission. 

2. 2. 2. 7  Summary  of  Gross  Channel  Parameters 

In  light  of  above  discussions,  the  gross  transmission  parameters,  important 
in  characterizing  the  effects  of  a  RTVL  medium  upon  broad  classes  of  signals,  are: 

1.  The  outage  (dropout  or  fadeout)  is  the  fraction  of  time  in  an  interval  in 
which  the  characteristics  of  the  useful  signal  parameters  do  riot  change  significantly. 
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The  signal  parameter  most  frequently  examined  is  the  envelope,  and  the  threshold 
is  determined  by  signal-relative-to-noise  considerations. 

2.  The  delay  spread  is  a  measure  of  the  duration  of  the  channel  impulse 
response,  or  of  the  maximum  delay  difference  between  the  first  and  last  slgnifi- 
cant  paths. 

3.  The  coherence  bandwidth  is  the  maximum  frequency  range  over  which  all 
frequency  components  of  a  signal  would  maintain  substantially  fixed  relative  ampli  ¬ 
tude  and  phase  relationships  and  fluctuate  practically  in  step  if  they  were  contained 
within  the  coherence  bandwidth.  The  channel  signal-distortion  effects  upon  link 
performance  would  thus  be  negligible.  The  bandwidth  is  inversely  proportional 

to  the  delay  spread. 

4.  The  diversity  (or  the  decorrelation)  bandwidth  is  the  minimum  frequency 
separation  between  two  Input  carriers  that  result  in  the  correlation  coefficient  of 
1/e  or  less  between  the  corresponding  envelopes.  The  diversity  bandwidth  is 
about  10  to  SO  times  the  coherence  bandwidth. 

5.  The  fading  rate  (or  Doppler  spread)  is  a  measure  of  the  nominal  width  of 
the  dispersion  in  frequency  experienced  by  each  individual  frequency  component  in 
the  transmitted  signal. 

6.  The  diversity  (or  decorrelation)  time  is  a  measure  of  the  time  separation 
that  must  exist  between  two  impulse  excitations  to  yield  a  correlation  coefficient 
of  1/e  or  less  between  the  output  envelopes.  The  decorrelation  time  is  inversely 
proportional  to  the  fading  bandwidth. 

7.  The  instantaneous  channel  function  h(t,  €)  which  is  discussed  earlier. 

8.  The  second-order  channel  functions  which  are  discussed  earlier. 

These  parameters  are  not  all  independent  because,  except  for  outage  probabil¬ 
ity,  they  are  determined  by  R^  (fi.  t)  a  td  SIt^.  f^)  of  Eqs.  (9)  and  (10).  Some  of 
the  interrelationships  are  already  discussed  earlier,  The  nominal  spread  of 
SlTj,  fj)  in  Tj  dimension  is  a  measure  of  the  multipath-delay  spread  of  the  multi¬ 
path  channel  or  of  the  delay  spread  oi  single  path.  The  nominal  spread  of  Sl-^.  f^) 
in  the  fj  dimension  is  a  measure  of  the  fpding  rate.  The  nominal  spread  of 
Rjjjj  (J>,  t)  in  the  O  dimension  is  r.  measure  of  (a)  the  low -distort  ion  transmission 
bandwidth,  when  the  delay  difference  between  two  paths  is  greater  than  the  delay 

u 

spreads  of  the  individual  path  (see  Figure  7)  ;  (b)  the  coherence  band- 

U 

width  of  paths  with  overlapping  delay  spreads  (see  Figure  8).  The  nominal 
spread  of  Rjjjj  (fi,  T)  in  the  t  dimension  is  a  measure  of  the  channel -decorrelation 
time. 


8.  l.omax,  J.B.  (1970)  ilF  propagation  dispersion,  in  Phase  and  Frequency 
Instabilities  in  Electromagnetic  Wave  Propagation.  K.  Davies,  editoi, 
pp.  497-510,  A<!j/iRD  Conference  Proceedings  No.  33,  Techvision  Services 
Slough,  England. 
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TIME - m» 

TIME  PROFILE 


Figure  7.  Scattering  Function  Showing  That  the  Delay  Difference  Between  Two 
Paths  Is  Greater  Than  the  Delay  Spreads  of  an  Individual  Path 


TIME  PROFILE  TIME - M* 

Figure  tt.  Scattering  Function  Showing  Paths  With  Overlapping  Delay 
Spreads 
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In  Section  2.4,  we  will  discuss  the  measurement  techniques  for  these  channel 
parameters.  In  the  following  section,  we  will  develop9  the  time -frequency  plane 
model  (or  "diversity  model")  used  In  signal  design  In  Section  2.  5,  from  the  scat¬ 
tering  function  S(Td,  f^). 

2. 2. 2. 8  Derivation  of  the  Time-Frequency  Plane  Model 
from  the  Canonical  Channel  Model 

In  this  section,  we  will  derive  the  time-frequency  plane  model  from  the  scat¬ 
tering  function  which  is  useful  in  signal  design,  we  will  also  demonstrate  the  var¬ 
ious  shapes  of  the  possible  delay-Doppler  occupancy  patterns. 

If  x(t)  is  the  Input  and  y(t)  Is  the  output,  then 


y(t)  *  //x(t  -  rd>  exp  J2»fd(t  -  td>  v  (rd,  fd)  drd  dfd  .  (13) 

where  V(rd,  fd)  =  J  h(vd,  t)  dt  is  called  the  Doppler-delay  spread  function. 

If  the  input  to  a  channel  is  confined  by  a  time  gate  to  the  time  interval 
0  <  t  <  T  and  the  output  spectrum  is  confined  by  a  bandpass  filter  to  the  frequency 
interval  *yW<u  <  y  W,  then  ln  of  the  actual  channel  with  Doppler-delay 
spread  function  V(rd,  fd)  one  may  use  a  channel  whose  Doppler-delay  spread  func¬ 
tion  v(rd,  fd)  has  the  singular  form 

rd’ s  l  l  v.,'  ('d-7){(Td-B)  1141 

m  ,  n 

where  6(  )  is  the  unit  impulse, 


//“" 

[M-V- 

)]  ainc  [T  ('a  ‘  f )] 

sine  [v 

'  (Td  • ® )] 

«'d'  V  dTd  ",j 

(15) 

and 


sine  x 


sin  *  x 
»  x 


9.  Bello,  P.A.  (1969)  Measurement  of  random  time-vai  ant  linear  channels, 
IEEE  Trans,  on  Info,  Th. ,  Vol.  IT-25,  No.  4,  July. 
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Examination  of  Eq.  (15)  shows  that  (V  }  are  essentially  two-dimensional 
sampled  versions  of  the  original  scattering  function,  the  sampling  taking  place 
with  "pulses"  of  width  of  the  order  of  1/T  in  the  fd  direction  and  1/W  in  the  Td 
direction. 

The  discrete  channel  model,  which  consists  of  a  finite  set  of  delays  and 

1  :  U|  <  T 

Doppler  shifts  is  shown  in  Figure  9,  where  rect(x)  -  7  .  The 

Oslxlssj 

output  y(t)  in  Eq.  (13)  can  therefore  be  written  as 


(Ui) 

It  is  known9  that  the  Doppler-delay  spread  function  V(Td,  f^)  and  the  scattering 
function  S(rd,  f^)  are  related  as 

V*(Td, 7d)  V<r“ >„)  -  S(Td,  fd)  (17) 


or 


s(S*t)/tw  *  f(,r  m  "  ;  3 

*  0  ,  for  m  r,  n*s 


(18) 


that  is,  the  scattering  function  varies  very  little  for  changes  in  rd  of  the  order  of 
l  /W  and  changes  in  fd  of  the  order  of  1  /T.  Thus  the  gains  of  the  discrete  point 
"scatte-ers"  become  uncorrelated,  and  the  strength  of  the  reflection  from  a  par¬ 
ticular  scatterer  becomes  proportional  to  the  amplitude  of  the  scattering  function 
at  the  same  value  of  delay  and  Doppler  shtft.  For  T  >•  than  the  multipath  delay 
spread  and  W  »  than  the  fading  bandwith,  the  number  of  V  coefficients  signif¬ 
icantly  different  from  zero  will  be  determined  by  how  muny  rectangles  of  limensinh 
1  •'>  W  can  be  fit  into  the  delay-Doppler  occupancy  pattern,  that  is  into  regions  of 
the  Tfj,  f^  plane  over  which  V(Tj,  f(|)  is  significantly  different  than  zero.  Such 
regions  arc  shown  in  Figure  10.  If  is  the  area  of  one  such  region,  then  the 
number  of  loefficients  of  significant  amplitude  can  be  expressed  as 


Ncoef  =  TW  SA 


(19) 
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2. 2.  2.  9  Noise  and  Interference 

The  levels  of  noise  and  interference  at  the  two  terminals  of  the  link  will  not 
be  the  same,  therefore  these  measurements  must  be  performed  at  both  stations. 

Receiver  front-end  noise  will  be  negligible  at  HF  with  respect  to  background 
noise  (atmospheric  and  galactic)  and  with  respect  to  man-made  interference. 
Table  2  provides  an  estimate  in  dB  above  KTB  of  the  atmospheric  and  galactic 
noise  level  at  20  MHz,  for  various  seasons  and  hours  of  the  day. 


Table  2.  Atmospheric  and  Galactic  Noise 
(dB  above  KTB) 


20  MHz 

Hours 

Atmo. 

Galactic 

Winter  0000-0400 

22 

22 

0400-0800 

24 

22 

0600-1200 

28 

22 

1200-1600 

42 

22 

1600-2000 

37 

22 

2000-2400 

30 

22 

Spring  0000-0400 

30 

22 

0400-0800 

25 

22 

0800-1200 

34 

22 

1200-1600 

45 

22 

1600-2000 

37 

22 

2000-2400 

39 

22 

Summer  0000-0400 

22 

22 

0400-0800 

25 

22 

0800-1200 

25 

22 

1200-1600 

36 

22 

1600-2000 

34 

22 

2000-2400 

30 

22 

Autumn  0000-0400 

22 

0400-0800 

22 

0800-1200 

36 

22 

1200-1600 

38 

2? 

1600-2000 

40 

22 

2000-2400 

36 

22 
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Interference  into  the  receivers  of  the  link  from  nearby  transmitters,  or  from 
co-channel  emissions,  either  at  close  range  or  at  a  distance,  as  well  as  man-made 
noise  will  be  the  predominant  factors  in  establishing  the  overall  signal-to-(noise 
+  interference)  ratio  at  a  given  carrier  of  the  HF  link.  As  it  will  be  discussed 
later  on  in  this  report,  the  adaptive  scheme  that  we  propose  will  use  a  waveform 
characterized  by  the  presence  of  numerous  spectral  lines  and  the  noise  3  interfer¬ 
ence  level  at  each  one  of  them  will  be  verified  in  order  to  exclude  the  most  inter¬ 
fered  ones. 

After  review  of  CCIR  Report  No.  65,  we  have  adopted  in  our  analysis  two 
values  for  the  level  of  noise  +  interference;  +30  dB  and  +*">  dB  above  KTB. 


2. 2.  2.10  Estimates  of  the  Path  Losses  and  of  the  Time/ 

Frequency  Spreads 

Ideally,  we  would  like  to  provide  reliable  estimates  of  all  the  parameters 
indicated  in  Figures  11,  12  and  13  for  each  path  of  interest  in  our  study.  They  are 
the  quantities  B^,  Bo,  Bm,  S(f,  v).  Actually  it  would  be  sufficient  to  simplify  the 
scattering  function  to  a  group  of  N  gaussoides,  with  B()  =  0,  B^  =  0,  and  to  reduce 
therefore  the  scattering  function  to  the  one  shown  in  cross-section  in  Figures  14 
and  15,  where  the  analytical  expression  of  S(5,  v)  becomes 


S(f,r) 


N 

£  P^tBI.,)'1  e\p 

i=  l 


(r  -  .v)2 


(20) 


In  this  formula,  the  parameter  N  represents  the  number  of  paths  in  the  struc¬ 
ture,  r j  and  L.  are  the  mean  delay  and  the  multipath  spread,  B  is  the  Doppler 
spread  of  the  path,  and  P(  denotes  the  relative  strength  of  the  i1*1  path. 

Further  simplification  can  be  achieved  by  representing  the  scattering  function 
S(5,  v)  as  a  single  gaussoid,  whose  amplitude  is  a  function  of  the  path  losses  and 
whose  widths  l.tot  and  B((>t  are  respectively  the  total  time  spread  and  the  total 
Doppler  spread. 

Because  of  the  time  variability  of  the  path,  we  also  require  the  knowledge  of 
the  statistical  properties  of  all  these  parameters,  so  that  <ve  can  compute  then- 
median  values  and  plot  the  curves  that  provide  ihe  percentage  of  the  time  (y-axis) 
during  which  path  losses,  the  multipath  spread  and  the  Doppler  spread  exceed  the 
value  of  the  abscissa  (x-axis).  Unfortunately  only  fragmentary  data  exist  that  are 
usable  to  this  end. 

First,  let's  review  the  case  of  a  one-hop  HF  path  m  midlatitude.  Davies^ 
gives  a  step-by-step  procedure  that  can  be  easily  followed  and  leads  to  a  reliable 
estimate  of  the  path  losses.  For  a  path  with  a  length  2.  5  megameters  working 

10.  Davies,  K.  (19tiS)  Ionospheric  Radio  Propagation,  NHS  Monograph  No.  80. 
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at  such  a  frequency  f  that  0.  85  MUF4000  <  f  <  MUF4000,  a  typical  value  of 
path  losses  (inclusive  of  absorption  losses  and  of  antenna  gains)  is  130  dB.  Few 
theoretical  formulations  and  even  fewer  experimental  data  are  available  on  the 
multipath  spread  and  Dopper  spread  of  single-hop  HF  paths  in  midlatitude.  Bailey*1 
investigated  the  HF  multipath  spread  phenomenon  and  its  dependence  upon  operating 
frequency  (specifically,  upon  its  ratio  to  the  MUF).  He  also  included  in  his  analysis 
the  dependence  of  multipath  spread  on  psth  length,  link  location,  local  time  t  3  well 
as  season.  This  analysis  shows  that  in  a  2.  5  megameter  path,  when  using  a  work¬ 
ing  frequency  of  0.  9-MUF4000,  the  time  spread  L  3  100  microsec.  Reliable  sta¬ 
tistics  on  the  time  spread  are  unavailable.  Concerning  now  the  Doppler  spread  of 
a  midlatitude  one-hop  HF  path,  we  can  estimate  a  median  value  of  0. 1  Hz,  For 

this  parameter  too,  a  reliable  statistic  is  an  unfulfilled  requirement. 

12  13 

Let's  review  now  the  case  of  a  transauroral  path.  Lomax,  Shaver  et  al, 

14 

and  Shepherd  and  Lomax  have  reported  on  experimental  measurements  of  high- 
latitude  HF  propagation  characteristics,  specifically  for  the  frequency  of  7.3G6  MHz. 
These  authors  found  that  winter  propagation  (inclusive  of  "normal"  and  "off-path" 
rays)  is  characterized  by  a  time  spread  of  approximately  1200  microseconds  and 
by  a  Doppler  spread  of  approximately  1  Hz.  Summer  propagation  shows  Doppler 
spread  up  to  20  Hz,  although  most  of  the  signal  energy  was  received  within  0.  3  Hz 
from  the  carrier.  Data  on  transauroral-path  losses  at  HF  are  scant,  although 
these  are  certainly  larger  than  in  midlatitude  paths.  In  our  study  we  have  assumed 
two  values  for  this  parameter,  150  dB  and  165  dB.  An  estimate  of  the  percentage 
of  time  during  which  these  values  will  be  exceeded  is,  however,  lacking. 

Table  3  summarizes  the  channel  properties  that  we  have  assumed  as  educated 
guesses  for  our  study  of  adaptive  HF  propagation  path  utilization. 


11.  Bailey,  D.  K.  (1959)  The  effect  of  multipath  distortion  on  the  choice  of 

operating  frequencies  for  HF’  communications  circuits,  IRE  Trans.  Antenna 
and  Propagation,  A  P-7,  398.  *  - 

12.  Lorrax,  J.B.  (1970)  HF  propagation  dispersion,  in  Phase  and  Frequency 

Instabilities  in  Electromagnetic  Wave  Propagation,  K.  Davies,  editor, 
pp.  497-510,  AGARD'  Conference  Proci  edings  No!  33,  Techvision  Services, 
Slough,  England. 

13.  Shaver,  H.N.,  Tupper,  B.C.,  and  Loma  J.B.  (1967)  Evaluation  of  a 

Gaussian  HF  channel  model,  IEEE  Tr  ,s.  on  Communications  Technology, 
Vol,  15,  No.  1,  pp.  79-88. - - 

14.  Shepherd,  R.A.  and  Lomax,  J.R.  (1967)  Frequency  spread  in  ionospheric 

radio  propagation,  IEEE  Trans,  on  Communications  Technology,  Vol.  15, 
No.  2,  pp.  268-2751 - 
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Table  3.  HF  Channel  Properties 


Midlatitude  Path 

Transauroral  Path 

Multipath  Spread 

100  psec 

1500  fisec 

Doppler  Spread 

0. 1  Hz 

10  Hz 

Path  Losses 

130  dB 

(150  dB 
\  165  dB 

Noise  +  Interference 

(30  dB 

(30  dB 

(above  KTB) 

\  45  dB 

l  45  dB 

2.3  PrUi  Sounding 

Path  sounding  has  the  scope  of  measuring  path  losses  at  an  adequate  number  of 
spot  frequencies  in  the  band  of  interest  (3  MHz  to  30  MHz)  and  of  measuring  at 
the  same  time  noise  and  interference  levels,  at  the  same  frequencies  and  at  both 
ends  of  the  link. 

Table  4  gives  the  parameters  of  the  proposed  sounding  scan.  The  master 
station  radiates  sequentially  1125  to  3750  carriers  to  cover  the  3 -MHz  t0  30-MHz 
band,  in  a  time  interval  100  to  160  seconds  (88  milliseconds  to  47  milliseconds  per 
carrier).  Of  the  two  numbers  given  above  for  each  sounding  parameter,  the  first 
applies  to  a  midlatitude  path,  the  second  to  a  transauroral  path.  The  scan  is  re¬ 
peated  every  5  to  8  minutes. 

Once  a  set  of  frequencies  has  been  chosen  for  communicating,  it  is  automat¬ 
ically  excluded  from  next  sounding  cycle.  However,  information  on  the  channel 
status  for  each  one  of  the  frequencies  thus  excluded  from  sounding  and  probing  is 
still  updated  once  every  5  to  8  minutes  by  measurements  oer formed  on  the  commun 
icat ion  waveform.  Frequency  switching  is  preceded  by  a  "tone"  of  notification  and 
takes  place  even  while  communications  go  on,  for  the  case  in  which  the  channel 
deteriorates  and  another  set  of  frequencies  is  found  more  suitable  for  carrying 
out  the  communications. 

The  block  diagram  in  Figure  16  has  been  worked  out  for  Link  Configuration  1 
(two-way  adaptive  link)  and  illustrates  the  various  functions  of  the  two  terminals  of 
the  link.  Here,  one-way  sounding  and  probing  is  achieved  from  the  master  station 
to  the  slave  station.  Acknowledgement  is  from  the  slave  to  the  master  station. 
Finally,  communication  is  a  two-way  exchange  between  the  stations. 
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Table  4.  Sounding/Probing  Scan  Parameters 


Midlatitude  Path 

Transaurora*  Path 

Band  covered 

3  MHz-30  MHz 

3  MHz-30  MHz 

Number  of  spot  frequencies 

1125 

3375 

Separation  between  two 
adjacent  spot  frequencies 

24  KHz 

8  KHz* 

Sounding  scan  time 

100  seconds 

160  seconds 

Bate  of  sounding  scan 
repetit’on 

one  every 

300  seconds 

one  every 

480  seconds 

Dwelling  time  per  spot 
frequency 

88  mlllisec 

47  milliaec 

Nominal  bandwidth  of 
sounding  receiver 

24  KHz 

8  KHz 

Width  of  sounding  pulse 

41.  5  microsec 

125  microsec 

Pulse  repetition  frequency 

100  pps 

100  pps 

Pulses  per  dwelling  time 

8  pulses 

4  pulses 

PRF  duty  cycle 

10*3 

1.25  10'2 

Noise  and  interference 
measurement's  integration 
time,  for  each  spot 
frequency 

53  mlllisec 

4.  millisec 

Overall  noise  and  inter¬ 
ference  measurement 
time 

60  seconds 

160  seconds 

$ 

This  value  is  chosen  because  £  KHz  is  the  bandwith  of  the  signal  waveftmm 
selected  for  the  transauroral  link.  The  path  cuhet .  »»t  bandwidth  is  only 
660  Hz. 


Midlatitude  link 


Step  1  -  100  seconds  devoted  to  sounding  operation. 

Step  2  -  SO  seconds  devoted  to  measurement  of  noise  and  interference  at 
both  terminals  of  the  link. 

Step  3-20  seconds  devoted  to  computations,  taking  into  account  the  need 

of  accumulating  at  a  single  terminal  (the  slave  station)  the  informa¬ 
tion  pertaining  to  noise  and  interference  at  both  terminals.  During 
this  step,  the  microprocessor  at  the  slave  station  selects  the  fre¬ 
quencies  and  designates  them  to  the  master  station. 


3  5 


3(» 


Figure  16.  Simplified  Block  Diagram  of  Link's  Terminals 


Step  4  -  100  seconds  devoted  to  channel  probing,  to  be  performed  only  at 
the  frequencies  designated  by  Step  3. 

Step  5-20  seconds  devoted  to  computations,  acknowledgement  and  informa¬ 
tion  exchange  between  the  two  terminals,  In  order  to  perform  the 
final  selection  of  frequencies  to  be  used  In  communications,  by  tak¬ 
ing  into  account  the  data  on  multipath  spread  and  Doppler  spread. 

Step  6  -  The  two  terminals  are  now  ready  to  initiate  communications.  The 
frequencies  finally  adopted  for  communications  are  excluded  from 
next  sounding/  probing  cycle  tone  every  300  seconds),  although  they 
continue  to  be  monitored  by  measurements  on  the  modulation  wave¬ 
form. 

Transauroral  l>nk 

Step  1  -  160  seconds 

Step  2  -  160  seconds 

Step  3-20  seconds 

Step  4  -  120  seconds 

Step  5  -  20  seconds 

Step  6  -  The  two  terminals  are  now  ready  to  initiate  communications. 

The  sounding/probing  cycle  is  repeated  every  8  minutes  (480  seconds). 

Communications  are  therefore  inhibited  only  in  the  first  300  (or  480  seconds) 
of  link  operation.  After  this  initial  adaptive  adjustments  of  the  link's  terminals, 
any  readjustment  is  performed  without  requiring  a  discontinuation  of  communica¬ 
tions. 

2.4  Channel  hubing 

The  importance  of  time-  md  frequency-dispersive  effects  in  HF  propagation 
has  been  amply  treated  in  the  literature  and  there  is  no  doubt  that  such  effects  are 
determining  factors  in  the  c  >nceptual  design  of  an  adaptive  system.  Channel  prob¬ 
ing  is  aimed  at  gathering  information  on  these  effects,  after  the  path  sounding  has 
determined  path  losses  and  no  ice  plus  interference  levels  at  th »  available  spectral 
lines,  and  has  identified  the  frequencies  promising  enough  to  be  worthy  of  the 
channel -probing  effort.  All  these  functions  are  slowly  varying  functions,  so  that 
one  sample  every  5  to  8  minutes  is  adequate. 

The  measurement  of  multipath  spread  and  of  the  Doppler  spread  can  be  achieved 
with  a  variety  of  methods,  either  based  on  the  direct  measurement  of  these  two 
quantities  or  on  indirect  measurements  such  as  the  ones  based  on  the  fact  that,  at 
a  given  frequency,  the  reciprocal  of  the  Doppler  spread  gives  the  e.  m.  wave-fading 
period  or  that  the  reciprocal  of  the  multipath  spread,  at  a  given  instant  m  time, 
gives  the  frequency  interval  within  which  car-iers  fade  coherently.  Because  the 


amount  of  time  required  to  process  the  Information  on  the  dispersive  properties 
of  each  channel  Is  not  trivial,  we  propose  to  perform  these  measurements  only  for 
those  frequencies  for  which  path  sounding  has  indicated  acceptable  path  losses  and 
affordable  noise  and  Interference  levels.  Therefore  channel  probing  has  to  follow, 
in  time,  the  sounding  operation,  as  was  illustrated  in  Section  2.3. 

In  the  following  sections  we  will  review  applicable  measurement  techniques. 

2.4.1  MEASUREMENT  TECHNIQUES  OF  TRANSMISSION 
PARAMETERS 

The  measure  of  communication  channels  is  important  in  digital  communications 

because  high-speed  digital  data  transmission  requires  considerable  knowledge  and 

3 

equalization  of  channel  characteristics.  Kallath  pointed  out  that  the  problem  of  the 

measurement  of  system  functions  (see  Section  2. 2. 2,  6)  of  random  time-variant 

channels  differs  from  the  corresponding  problem  for  time-variant  channels  in  that, 

even  in  the  absence  of  noise,  the  random  system  function  may  be  unmeasurable. 

Kallath  introduced  a  channel  parameter  called  a  "spread  factor"  as  the  measurability 

criterion.  This  parameter  is  the  product  of  Bmax,  the  maximum  rate  of  variation 

of  the  system  in  Hz,  and  Lmax,  the  maximum  multipath  spread  of  the  channel  in 

seconds.  According  to  Kallath,  the  system  ftmctlons  of  a  linear  channel  cannot  be 

measured  if  the  "rectangular  spread  factor*'  of  the  channel,  S„  (=  B„„„  L  „)  >  1 

It  max  max 

and  if  no  further  information  than  13  L _ Is  known  about  the  channel.  How- 

max*  max 

ever,  he  was  careful  to  point  out  that  additional  channel  knowledge  would  generally 

allow  exact  channel  measurement  even  though  B  „  L _ >  1.  Based  upon  the 

max  max 

discrete  representations  of  the  channel  (Figure  9),  corresponding  to  input -time  and 

o 

output -bandwidth  constraints,  Bello  proposed  the  region  of  the  r.on-zero  delay- 
Doppler  occupancy  pattern  (Figure  10)  being  less  titan  unity  as  the  less  stringent 
measurability  criterion.  This  new  channel  parameter  is  called  the  area  spread 
factor  of  the  channel  S^.  Another  method  of  evaluating  the  'spread  factor'  of  the 
channel  is  to  define  the  spread  factor  r.a  the  ratio  of  the  bandwidth  of  the  fast  fluc¬ 
tuations  and  the  correlation  bandwidth  of  the  channel  being  measured.  Fortunately 
except  for  certain  esoteric  radio  channels  (for  example,  the  orbital  radio  channels), 
the  spread  factor  is  less  than  unity  and  thus  the  measurement  techniques  described 
below  ate  generally  quite  useful. 

These  measurement  techniques  for  a  random  dispersive  channel  are  analyzed 
in  three  levels  of  increasing  complexity: 

1.  The  measurement  of  multipath  spread  and  Dopper  spread,  and  Doppler  shift 
and  spectral  skewness. 

2.  the  measurement  of  second-order  channel  functions: 

3.  the  measurement  of  instantaneous  channel  functions. 

For  the  parameters  in  (l),  measurement  techniques  used  are  based  upon 
differentiation,  level-crossing  and  correlation.  For  (2),  the  techniques  used  are 
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correlation  technique,  multitone  technique,  pulse-pair  technique  and  the  chirp 
technique.  For  the  measurement  o f  the  Instantaneous  valves  of  the  channel  func¬ 
tions,  the  cross-correlation,  multitone  and  pulse-patr  techniques  are  used. 

2.4.2  MEASUREMENT  OF  DOPPLER  SPECTRUM 
PARAMETERS 

In  this  section,  we  will  summarize  the  techniques  of  the  real-tir  measure¬ 
ment  of  the  centroid,  rms  bandwidth  and  skewness  of  the  Doppler  spectrum  of  the 
complex  envelope  of  the  narrowband  process.  Due  to  the  Doppler  spreading  char¬ 
acteristics  of  the  channel,  the  received  process  is  narrowband  when  a  CW  tone  is 
transmitted.  The  importance  of  the  rms  bandwidth  (or  second  central  moment  of 
the  spectrum)  in  signal  design  and  in  probability  of  error  computations  is  well 
known  in  the  literature.  In  Appendix  B,  we  have  outlined  the  effect  of  centroid 
(first  moment)  and  the  skewness  of  the  spectrum  on  the  probability  of  error  com¬ 
putations  and  the  equivalent  signal-to-noise  ratio  (S/N).  This  S/N  parameter  is 
important  in  evaluating  the  number  of  frequencies  required  in  combatting  the  fading 
and  inter-symbol  interference  (see  Section  2.  5).  In  that  section,  the  probability  of 
error  expression  only  accounts  for  the  additive  noise  in  the  channel  and  neglects 
the  fading  of  the  channel.  However,  we  propose  lhat  by  using  the  equivalent  signal- 
to-notse  ratio  in  the  probablllty-of-error  expression,  the  fading  effect  of  the  chan¬ 
nel  can  also  be  accounted  for.  Gupta15  has  recently  developed  techniques  to  com¬ 
pute  the  spectral  moments  of  the  complex  envelope  of  the  narrowband  process. 

The  main  results  are  summarized  here: 

"General  expressions  of  the  spectral  moments  of  any  order  (specitlcally  the 
first  four)  are  derived  by  the  author  in  the  time  domain  and  their  derivations  do  not 
require  the  computation  of  the  power  spectrum.  Ir.  order  to  simplify  the  signal 
processing  of  the  spectral  momei.’  estimators,  the  author  proposed  a  general  class 
of  these  estimators  of  an  undistoited  complex  envelope  which  are  obtained  from 
one  or  both  distorted  quadrature  components  of  the  complex  process.  The  com¬ 
ponents  are  distorted  by  means  of  mentoryless  nonlinear  devices  (MLNLD),  such 
as  the  hard-clipper.  By  making  use  of  the  previous  results  of  Bussgang  (cross - 
correlation  theorem)  and  Van  Vleck  (arcsln  law),  interesting  results  are  obtained. 
For  example,  it  is  shown  that  odd  (even)  spectral  moments  of  an  undistorted  com¬ 
plex  process  can  be  obtamed  when  both  quadrature  components  are  distorted  by 
even  (odd)  MLNLD  and  that  even  (odd)  moments  cannot  be  obtained  if  both  compo¬ 
nents  are  hard-clipped  (even-linearly  rectified).  However,  the  spectral  moments 
of  any  order  can  be  computed  in  the  time  domain  when  only  one  of  the  two  compo¬ 
nents  is  distorted  by  arbitrary  MLNuD." 


15.  Gupta,  A.K.  (1979)  Dissertation  abstract,  IEEE  Trans.  Info.  Theory,  p.  760  , 
November. 
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Mathematically,  if  S(f)  is  the  power  spectrum  of  a  received  carrier,  the  cen¬ 
tral  spectral  moments  are  defined  as 


mn 


fit  -  S(f)  df 
f  S(f>  df 


(21) 


and  the  simple  spectral  moments  MR  are  obtained  by  replacing  (f  -  Pj)n  by  f”.  It  is 
simple  to  obtain  the  relationship  between  pn  and  mfl.  is  called  the  centroid, 

m„  and  m,  are  second  and  third  central  spectral  moments  respectively.  Skewness 

^  J  15 

is  defined  by  Gupta  as 


s 


(22) 


If  g(t)  is  the  complex  envelope  of  the  received  process  and  is  given  by  g(t)  *  o(t)  + 
jfl(t),  then  for  the  case  of  undistorted  components 


the  rms  bandwidth  or  Doppler  spread  (D)  is  given  by  twice  vAng  .  a  and  8  are  o(t) 
and  flit)  respectively,  denotes  differentiation  and  "{  )"  denotes  ensemble 
averaging. 

The  inphase  and  quadrature  components  a(t)  and  (3(t)  can  be  determined  by 
multiplying  the  received  carrier  by  both  a  local  carrier  and  a  90°  shifted  local 
carrier  at  the  same  frequency  as  the  received  carrier  (or  as  near  to  the  same 
frequency  as  possible)  and  then  extracting  the  low-frequency  components.  Strictly 
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speaking  JJ(=2Vm,)  and  skewness  (s^)  are  independent  of  centroid,  and  thus  pre¬ 
cise  knowledge  of  the  received  carrier  frequency  is  not  necessary.  However,  as 
the  local  carrier  frequency  departs  from  the  received  carrier  frequency,  the 
extracted  o( t)  and  /5(t)  Increase  in  bandwidth,  necessitating  larger  bandwidth  filters 
and  passing  more  noise.  Thus,  from  the  point  of  view  of  maximizing  signal-to- 
noise  ratio,  it  is  desirable  to  keep  the  local  carrier  frequency  as  near  as  possible 
to  the  received  signal  frequency. 

We  consider  now  a  simpler  technique  for  the  measurement  of  Doppler  spread 
that  uses  only  the  envelope  or,  more  generally,  any  well-behaved  nonlinear  function 
of  the  envelope  of  the  received  carrier.  For  the  technique  proposed  to  be  strictly 
correct,  It  Is  necessary  to  assume  that  the  transmission  of  a  carrier  results  in 
the  reception  of  a  narrowband  Gaussian  process.  However,  slight  departures  from 
Gauss ianity  should  not  affect  the  measured  parameter  significantly. 

7 

It  is  demonstrated  that  If  e(t)  Is  some  nonlinear  fuctlon  (K)  of  the  envelope 
o  2 

a  +  0  of  the  received  carrier,  the  rms  Dcnpler  spread  is  then  given  by 


d  .  i.  /  «*<*»;*> 

9a  V  «e<t»2> 


(27) 


where  a  Is  a  constant  dependent  upon  the  nonlinear  device  and  is  r.iven  by 


a  = 


’r  k2(r)  dr 


(28) 


In  the  case  of  a  linear  envelope  detector  a  =  1/V2  and  for  a  square-law  detector 

o-l.  The  formation  of  the  derivatives  of  the  envelope  from  the  sampled  data 

1  fi 

requires  some  care.  In  Bello,  some  attention  has  been  given  to  this  problem. 

Fading  rates  have  also  been  determined  by  measuring  the  average  number  of 
times/unit  time  the  envelope  of  the  received  carrier  crossed  a  sp.  cified  level. 
Rice  (see  ref.  in  Bello  }  has  shown  that  for  a  narrowband  Gaussian  process  the 
number  of  times  (n)  the  envelope  crosses  a  level  R  is  simply  related  to  a  rms 
bandwidth  measure.  If  the  level  is  the  rms  of  the  envelope,  then 


16.  Bello,  P.  A.  (1965)  On  the  rms  bandwidth  of  non-ltnearly  envelope  detected 
narrowband  Gaussian  noise,  IEEE  Trans,  on  Info.  Theory,  pp.  236-239. 
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Another  approach  is  to  compute  the  zero  crossings  of  the  inphase  (or  quadrature) 
component  of  the  narrowband  Gaussian  process.  If  m  is  the  number  of  zero  cross¬ 
ings  then 

m  =  D  .  (30) 

Various  methods  of  estimating  centroid,  skewness,  and  Doppler  spread  are  dls- 
15 

cussed  by  Gupta. 

2.4.3  MULTIPATH  SPREAD  MEASUREMENTS 

The  multipath  spread  parameter  is  a  measure  of  the  dispersion  in  path  delays 
suffered  by  a  process  propagated  through  a  random  channel  (see  Figures  3,  7  and 
8).  If  Q(?)  is  the  delay  power  spectrum  which  describes  the  distribution  of  power 
in  the  various  path  delays,  the  rms  multipath  spread  is  given  by  (similar  to 
Eqs.  (21)  and  (26)) 


is  the  centroid  orQ(€)  (similar  topj  in  Eq.  (21)).  Although  the  following  analysis 
does  not  differentiate  the  discrete  scattering  function  (Figure  7)  and  continuous 
scattering  function  (Figure  8),  the  future  work  will  develop  new  methods  of  defin¬ 
ing  multipath  spread  in  two  different  situations. 

From  a  strictly  mathematical  point  of  view,  the  multipath-spread  measure¬ 
ment  problem  is  entirely  analogous  (dual)  to  the  Doppler-spread  measurement 
problem.  The  function  Q(?)  is  the  "power  spectrum"  of  the  instantaneous  channel 
transfer  function  G(f).  Thus,  if  a  very  narrow  pulse  of  spectrum  H(f)  is  trans¬ 
mitted,  W(f)  the  complex  amplitude  spectrum  of  the  received  transient  is  found  to  be 

W(f)  =  H(f)  G(f)  .  (33) 

If  H(f)  is  chosen  to  be  essentially  constant  over  a  band  of  frequencies  sufficient¬ 
ly  large  compared  to  the  correlation  bandwidth,  then  sufficient  values  of  G(f)  can 
be  determined  from  W(f)  to  perform  frequency-domain  averaging.  If  similar  to 


g(t),  one  defines  G(f)  as  G(f)  =  o(f)  +  jf3(f),  then  one  obtains15  from  Eqs.  (24)  and 
(26)  that 


where  "{  }"  denotes  frequency-domain  averaging  similar  to  the  relationships  in 
Eqs.  (23)  through  (25).  The  authors  of  this  report  point  out  that  M-  can  also  be 

a 

simplified  to 


(35) 


Similar  to  Eq.  (27),  one  can  also  obtain  M,,  from  the  envelope  as 


2  2 

where  E(f)  is  some  nonlinear  function  K(o  (f)  f  Jl  (f))  of  the  frequency  squared 
envelope. 

Similar  to  Eq.s.  (29)  and  (30),  multipath  spread  can  also  be  obtained,  from  (by 
using  time-frequency  duality)  the  average  number  of  times  per  hertz  the  transfer 
function  envelope  G(f)  crosses  its  rms  value  and  from  the  number  of  times  per 
hertz  the  real  or  imaginary  part  of  the  channel-transfer  function  crosses  zero. 

2.  4.  4  SIMULTANEOUS  MEASUREMENTS  OF  DOPPLER 
SPREAD  AND  MULTIPATH  SPREAD 

In  earlier  sections,  techniques  for  the  instantaneous  measurement  of  Doppler 
and  multipath  spread,  which  require  either  the  extraction  of  complex  envelopes  or 
envelopes  of  the  received  carriers.  In  this  section,  we  will  simultaneously  meas¬ 
ure  both  the  parameters  from  the  envelopes  of  the  received  carriers.  FM  and  SSB 
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17 

(single  sideband)  techniques  have  been  discussed  in  Bello.  In  this  section  we 
will  discuss  the  SSB  technique  only. 

In  the  SSB  technique  two  carriers  separated  by  F  Hz  are  transmitted.  Letting 
the  received  detected  envelopes  be  represented  as  Dt  (f  +  F)  and  Et(f),  the  normal¬ 
ized  difference  in  Eq.  (36)  is 


dEt(f)  Ef(f  +  F)  -  E(f) 
~nr~  =  F 


(37) 


where  subscript 't'  denotes  that  E(f)  is  also  function  of  time.  Assuming  stationary 
scatter  statistics,  time  and  frequency  averaging  can  be  interchanged.  Thus  in 
Eq.  (36) 


Figure  17  shows  the  block  diagram  implementation  of  Eqs.  (27)  and  (36)  with 
Eq.  (3 B)  using  the  SSB  technique.  Note  Mg  and  O  are  simultaneously  measured. 


Figure  17.  Simultaneous  Measurement  of  Doppler  Spread  and  Multipath  Spread 
Using  Envelopes  Only  by  the  SSB  Technique 


17.  Bello,  P.  A,  (1965)  Some  techniques  for  the  instantaneous  real-time  measure¬ 
ment  of  multipath  and  Doppler  spreads,  IEEE  Trans,  on  Comm.  Tech., 

Vol.  13,  No.  3,  pp.  155-192.  '  “  ' 


44 


2.4.5  MULTIPATH  MEASUREMENTS 


2.4.  5.1  Incoherent  Pulse  Probe 

A  suitable  signal  for  entir  .ating  the  delay  power  spectrum  Q(?)  is  a  short 
pulse  of  width  A,  which  is  .'arrow  compared  with  the  characteristic  variations  in 
Q(5).  Consider  the  signal, 


z(t)  =  p(t) 


p(t)  = 


0  |t|  >  A/2 

1  |t|  <  A/2 


A  <  mode  width  (for  example,  5  psec) 


(39) 


An  estimate  can  be  formed  by  square-law  detecting  the  observed  process  and  sub¬ 
tracting  off  a  noise  bias.  Thus,  for  the  response  to  one  of  the  pulses  in  the  train 
we  have 

Q(t)  *  |  r(t)| 2  -  |n(t)|2 

5  kA  J  g(£)  p(t  -  C)  dl  ¥  n(t) | 2  -  on  (40) 

The  properties  of  this  estimator  can  be  listed  as  follows: 

Mean 

Q(t)  =  A  P(t  -  €j)  g*(C2)  p(t  -  ?2)  d€j  dC2  +  | n(t)| 2  -  on 

*  A  J Q(C)|  p(t  -  ?)|2  d? 

a  EQ(t)  (41) 

where 

E  is  the  equivalent  pulse  energy,  that  is, 

E  •  A  J  !pvt>l2  dt  *  AA 

Thus,  Q(t)  will  be  a  filtered  version  of  the  true  function  if  enough  averaging  is 
carried  out.  The  resolution  is  deterr':ned  by  the  signal  properties,  so  that  the 


45 


pulse  width  must  be  narrower  than  the  mode  widths  of  Q(C)  to  give  an  adequate 
estimate. 

Variance 

Since  r(t)  is  a  complex  Gaussian  process, 

!  r(t)| 4  =  2[|  r(c> |  2] 2 
Thus, 


Var  (Q(t) }  =  ( j  r(t)| 2  -oj2  -  [Q]2 

■  |  r(t>  | 4  -2on  |  r(t)| 2  +  a2  -  (<§) 

2(E<5(t)  +  an) 2  -  2on  [E§(t)  +  aj  +  a2  -  (  E$(t)] 2 

*  (  E$(t)l 2  +  2onE$(t)  +  o2 

The  ratio  of  standard  deviation  to  mean  is 

i.d,  (Q(O)  s  1  EQ(€)  f  °n] 
mean  ((3(C))  eQ(C) 

. ,  ._f&_ 

eQ(0 


(42) 


(43) 


It  is  a  simple  matter  to  show  that  when  N  independent  "snapshots"  of  the 
channel  response  are  processed  in  this  way  and  averaged  by  overlapping,  the  ratio 
of  standard  deviation  to  mean  is  reduced  by  l/'^Ngq.  that  is, 


a.  dev. 
mean 
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For  a  required  ratio  of  0.  2  this  implies 


>  25 


Note  that  even  in  the  absence  of  noise  there  is  a  basic  limitation  to  the  attain¬ 
able  accuracy  for  a  given  measurement  time.  The  maximum  number  of  independent 
pulse  responses  which  can  be  obtained  in  1  minute  is 


N 


eq 


*  60 


Btot  -  Doppler  spread 


(44) 


The  ratio  of  standard  deviation  to  mean  for  this  zero-noise  case  is  found  to  be: 


a.  dev.  (Q(€))  _  1 

mean  (Q(5))  n/N_  _ 

vv| 

«  0.  13 


(45) 


In  the  next  few  sections,  which  include  additive-noise  calculations,  a  performance 
index  p  will  be  used,  where 


P  * 


EC$(€)  / 


(46) 


Examples 

Consider  a  multimodal  form  for  Q(?)  with  M  equal  energy  modes  of  width  Lq 
(see  Figure  13).  Then 

<5($)  *  (for  each  mode)  (47) 
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Since 


(48) 


Because  the  noise  level  varies  with  frequency  and  several  different  multipath 
structures  are  expected,  it  will  be  necessary  to  consider  several  sets  of  param¬ 
eters. 

Some  allowance  must  be  made  for  accurate  measurement  of  the  low-level 
detail  of$(?),  so  that  this  single  pulse  scheme  appears  unsuitable  because  of  low 
signal -to-noise  properties. 

2, 4.  5. 2  Coherent  Processing  Using  a  Coded  Sequence 

An  alternative  estimate  scheme  follows  if  coherent  processing  can  be  carried 
out  at  the  receiver.  Let  z(t)  be  a  pseudo-random  sequence  of  period  T  which  is 
used  to  phase  modulate  the  carrier  (bit  length  -  A). 

When  the  sequence  is  periodic  it  can  be  assumed  to  have  the  following  useful 
property 


T 

J  z(t)  z*(t  -  f)  dt  =  6(f) 
0 


with 


6(f) 


6(o) 
<M  o) 

s^N 


T  =  0 


2?  A 


(49) 


(!>0) 
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where 


N  =  sequence  length 
=  T/A 
=  2n  -  1 


(51) 


and  $(o)  is  the  energy  in  one  period 

However,  we  are  interested  in  transmitting  only  a  single  period  before  switch¬ 
ing  to  a  different  channel.  Thus,  for  a  particular  channel 


0 

z(t) 


0  •>  t  >  T 

otherwise 


(52) 


The  autocorrelation  for  this  aperiodic  zQ(t)  sequence  has  been  shown  to  have 
properties  similar  to  those  for  a  periodic  sequence,  that  is, 


6  It)  •  f  z  (t)  z  (t  -  t)  dt 

J  o  O 


-oe 

♦0(°) 

T  =  0 

*o<°> 

IV 

t> 

«/N 

We  now  consider  an  estimate  of  the  form 


Q(t) 


/ 

-00 


2 

r(t)  z  *  (t  -  7)  dt 

o 


b 


b  is  a  bias  term  which  will  be  discussed  subsequently.  Then 


(54) 
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00 

Q(r)  =  A  J  g(5)z0(t  -5)  g*(T)  zo*(t  -  t)  d?  dt2 
-00 


+ 


n(t)zQ*(t 


2 


t)  dt 


b 


=  A  | /g«H0(r  -  f)  d?  j 
00 

*  A  J  Q<€ )  k0<T  -  5 > | 2  d€  (55) 

-00 


b  ■ 


00 


J  n(t)  zo*(t  -  t) 


dt 


-oo 


o, 

iv 


I  /  iv(t-T>r 


dt 


-oo 


Aa^0(o) 


(56) 


where  b  has  been  chosen  so  that  the  second  term  in  the  expression  above  is  re¬ 
moved.  Using  the  relationf  hips  of  Eqs.  (55)  and  (53)  it  can  be  seen  that 


§<t)  =  A  fQ(5)U  (t  -  |)|2  df 

w  U 

*  A  J  Q(? ) |  p(S  -  v)|2  df  +  CQ 

*  AA«2(0)Q(5)  +  Oo  (57) 

where  C  is  a  bias  term  d".e  to  the  sidelobe  structure  of  Id  (t)|2,  and  |p(T)|2  rep- 
resents  |d(TM  in  the  vicimiy  of  the  origin  (that  is,  for  |t|  <  A). 
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§{!•)  as  it  is  defined  above  is  the  convolution  of  Q(?)  with  a  pulse  of  unit  area, 
so  that  this  notation  is  consistent  with  previous  analysis.  The  bias  term  CQ  can  be 
bounded  using  Eq.  (53). 


CQ  *  A  J  Q(f)  ^ 


M2J  0) 


o'u/  r 
f  Q«) 


A*J(0) 


For  a  true  function  Q(?)  with  modes  of  height. 


Q«)  - 


the  ratio  of  the  average  measured  mode  height  to  the  bias  CQ  will  be  of  the  oider, 


g<T)  aa^(Q)/mlo 
"'o  A$2(0)/N 


=  (NA/Ml.J 


The  worst  bias  effects  will  occur  when  A(t)  is  spread,  so  that  the  modes  are  of 
low  magnitude. 

Finally,  it  should  be  observed  that  although  it  is  them  o'ically  possible  to  re¬ 
move  this  bias  by  subtraction,  since  it  appears  to  depend  only  on  known  parameteis, 
in  practice  the  bias  will  be  a  function  of  r  because  of  variations  in  the  sidelobe 
structure  of 

Variance 

4  I  2 

VarlQU)]  =  Jr(t)zo*(t  -  ?)  dt  -  2b  |  fr(t)zo*(t  -  ?)  dt  +  b2  -  (Q) 2  (Gl) 
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Denote 


E  =  A 


f  U(0|2 


dt 


«  A*  JO) 

=  AT  (66) 

Note  that  in  computing,  the  CW  power  output  of  the  transmitter  must  be  used, 
rather  the  peak  power, 

2.4.  5.3  IMPLEMENTATION 

The  integral  of  Eq.  (54)  can  be  implemented  by  multiplying  the  original  se¬ 
quence,  which  can  be  generated  using  phase-reversal  keying  thus  implying  tnat 
z(t)  is  real,  against  the  inphase  and  quadrature  components  of  the  received  wave¬ 
form.  This  is  followed  by  a  square-law  detector  (or  possibly  linear  detection) 
and  an  integrate  and  dump  procedure  using  a  low-pass  filter  with  a  time  constant 
greater  than  T.  To  implement  the  scheme  for  a  set  of  >  values  requires  parallel 
processing  and  possibly  the  use  of  a  tapped  delay  line.  The  number  or  taps  depends 
on  the  product  of  the  signal  banawidth  and  the  range  of  delay  (r)  to  be  observed. 
Figures  18  to  22  indicate  the  basic  structure  of  the  receiver  and  of  the  transmitter. 
Not  included  is  the  logic  required  to  switch  through  the  set  of  frequencies  to  be 
measured. 

2.4.  5.4  Doppler  Spectrum  Measurements 

The  Doppler-spreading  characteristics  of  the  channel  are  best  observed  using 
a  CW  source.  However,  the  fading  rate  is  relatively  low,  so  that  efficient  Doppler 
measurements  on  a  set  of  chtnnels  requires  the  use  of  a  sequential  sampling 
approach.  Earlier  in  this  report,  it  was  noted  that  individual  channels  had  to  be 
probed  periodically  at  intervals  of  Tq  tTQ  >  l/Btot>  to  give  useful  estimates.  By 
probing  at  a  faster  rate  than  required  (T  <  l/BtQl),  information  is  obtained  that  is 
redundant  in  terms  of  the  delay  power  spectrum  estimate  (or  signal  level)  but  which 
may  be  useful  for  Doppler -sped rum  estimation. 

Consider  the  sequence  of  outputs  (f(nTQ)}  obtained  when  Eq.  (54)  is  tmpie- 
mtnted.  For  a  particular  r, 

f(nTQ)  =  QnT  (r)  (estimate  at  nTQ) 
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167) 


REFERENCE 


Figure  18.  Transmitter  System 


Figure  19.  Receiver  Overall  Block  Diagram 
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Figure  20.  AFC  Subsystem 
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REGISTER 


Figure  21.  Demodulators  I  and  Q 


ADVANCE 

COMMAND 


RETARD 
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TO  ADDITIONAL 
'SHIFT  REGISTER 


Figure  22.  PRN  Generator  and  Shift  Registers 
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Then  it  can  be  shown  that  if  the  sidelobe  interference  is  neglected, 

f(nT0)  *  |g(nTo,r)|2  (68) 

Thus,  any  spectral  information  obtained  from  this  sequence  will  relate  to  the 
shape  of  the  scattering  function  for  this  particular  delay  t. 

Because  of  the  inherently  higher  signal-to-noise  ratio,  it  appears  that  the 
latter  sequence  (using  simple  on-off  pulses)  offers  the  best  alternative  for  spec¬ 
trum  analysis.  In  addition,  the  total  scan  time  Tq  required  to  probe  all  channels 
once  is  significantly  less  than  for  the  wideband  coded -sequence  measurement,  thus 
providing  a  higher  sampling  rate  and  less  danger  of  aliasing  effects. 

It  should  be  noted  that  the  phase  Instability  of  the  transmitter  oscillator  over 
time  intervals  of  order  Tq  necessitates  the  use  of  an  incoherent  Doppler  measure¬ 
ment,  If  phase  coherence  of  the  oscillator  could  be  assumed,  then  in-phase  and 
quadrature  components  generated  at  each  Tq  seconds  would  be  used  to  get  the 
spectra1  properties  directly. 

2. 4.  6  MEASUREMENT  OF  INSTANTANEOUS  IMPULSE 
RESPONSE 

In  this  section,  we  discuss  the  measurement  of  the  instantaneous  values  of  the 

17 

impulse  response  of  a  continuous  multipath  channel  when  a  periodic  train  of 
identically  shaped  pulses  are  transmitted.  The  time  duration  of  the  typical  pulse 
is  chosen  small  enough  to  provide  the  required  resolution  to  characterize  the 
channel  for  signals  of  a  specified  bandwidth  and  should  also  be  less  than  the 
reciprocal  of  the  Doppler  spread.  The  pulse  duration  should  also  be  greater  than 
the  sum  of  the  channel-impulse  response  and  the  transmitted  pulse,  so  as  to  allow 
a  typical  pulse  response  to  become  very  small  before  a  new  pulse  arrives. 

If  the  complex  envelope  of  the  impulse  response  of  the  channel  is  denoted  by 
g(t,  ?)  and  the  additive  noise  by  n(t),  then  the  complex  envelope  of  the  received 
process  is  given  by 

w(t)  =  L  f  p(t  -  kT  -  5)  g(t,  5)  d?  +  n(t)  (69) 

k  J 

where  p(t)  denotes  the  complex  envelope  of  the  typical  transmitted  pulse.  To  charac¬ 
terize  a  channel  for  use  with  a  signal  of  bandwidth  VV,  it  is  only  necessary  to  meas¬ 
ure  a  smeared  g(t,  ?),  g(t,  S)  corresponding  to  the  channel  response  for  a  pulse  p(t) 
whose  spectrum  is  flat  over  the  bandwidth  W,  that  is, 

g(t,  n)  =  j p(n  -  f)  g(t,  5)  d?  .  (70) 
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If  the  received  process  w(t)  in  Eq.  (69)  is  sampled  with  the  periodic  impulse 
train  of  complex  envelope  T  =  £  6(t  -  n  -  k  -  T),  then  sampled  w(t)  is  given  by 

Tw(t)£6(t  -  n  -  kT)  =  T£6(t  -  n  -  kT)  g(t,  n)  +  TE6(t  -  rj  ••  kT)  rj(t)  .  (71) 

Thus,  the  signal  component  in  the  complex  sampled  version  of  the  received  process 
is  a  sampled  version  of  the  desired  impulse  response.  Complex  low -pass  filtering 

o 

will  then  recover  the  desired  impulse  response.  It  is  shown  in  Bello  that  for  a 
small  spread  factor  of  the  channel,  the  impulse  response  can  be  measured  with  a 
fairly  low  signal-to-noise  ratio  of  the  received  r.  f.  signal. 

Another  method  of  measuring  the  impulse  response  of  the  channel  is  the  well- 
known  cross -correlation  technique  which  involves  a  computation  of  the  short-time 
cross-correlation  function  between  the  output  of  the  channel  and  the  input  of  the 
channel.  The  averaging  time  is  adjusted  to  be  long  compared  to  the  time  constant 
of  the  transmitted  probing  signal  but  short  compared  to  the  fading-time  constant  of 
the  channel.  To  form  the  short-time  cross-correlation  function,  the  received 
signal  is  multiplied  by  a  shifted  replica  of  the  conjugate  of  the  transmitted  signal. 
The  desired  channel  Information  is  then  extracted  by  filtering  the  product,  using  a 
bandwidth  equal  to  the  fading  bandwidth. 

If  the  multipath  structure  of  the  channel  is  either  discrete,  or  for  nil  practical 
purposes,  may  be  regarded  as  discrete  for  the  class  of  input-band  limited  signals 
for  which  the  channel  is  to  be  characterized,  the  impulse  response  of  such  a  chan¬ 
nel  can  also  be  measured  by  the  above-mentioned  techniques,  If  a  train  of  short 
pulses,  short  enough  to  resolve  the  individual  contributions  of  the  discrete  paths, 
is  transmitted,  the  receiver  matches  ihe  received  process  by  a  shifted  trans¬ 
mitted  pulse  train.  An  integrator  notes  at  which  values  of  shift  the  maximum  out¬ 
puts  occur,  and  thus  determines  the  set  of  discrete  delays.  Once  these  are  deter¬ 
mined,  the  separate  paths  can  be  extracted.  Similar  to  the  continuous  delay  case, 
the  disci'o'e  paths  enn  also  be  estimated  by  the  cross-correlation  technique. 

2.4.  7  MULTITONE  MEASUREMENT  OF  THE  TIME- 
VARIANT  TRANSFER  FUNCTION 

The  transfer  function  of  the  channel  can  be  constructed  from  frequency -domain 
samples.  A  probing  signal  consisting  of  tones  spaced  DHz  apart  are  transmitted 
and  are  later  extracted  from  the  corresponding  received  fading  tones.  The  tone 
spacing  should  be  greater  than  the  Doppler  spread  to  obtain  spectrally  non-over¬ 
lapping  adjacent  tones.  The  tone  spacing  should  also  be  less  than  the  reciprocal 
of  the  multipath  spread. 
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2.4.8  MEASUREMENT  OF  CHANNEL-CORRELATION 
FUNCTIONS 

In  Section  2.  2.  2,  we  outlined  three  channel-correlation  functions:  tap-gain 
correlation  function  R(S,  t);  scattering  function  S(t^,  f^);  and  time -frequency  corre¬ 
lation  function  R(S7,t).  If  the  channel  is  represented  as  a  WSSUS  channel,  then  the 
purpose  of  most  channel-probing  techniques  is  to  measure  one  or  more  of  these 
three  channel-correlation  functions  and  associated  parameters.  These  probing 
techniques  are  designed  to  measure  the  average  fading  dispersive  characteristics 
of  a  link.  In  this  section,  we  discuss  four  channel-probing  techniques,  namely  the 

correlation  technique,  the  multitone  technique,  the  chirp  technique  and  the  pulse  - 

7 

pair  technique.  It  has  been  shown  in  Bello  that  all  these  techniques  yield  essen¬ 
tially  the  same  performance  when  the  measurement  error  is  minimized  for  a 
given  transmitted  energy.  The  ratio  of  the  square  of  the  maximum  value  of  the 
chaunel-correlatlor.  function  to  the  mean-square  error  in  the  measurement  of  that 
channel-correlation  function  is  given  for  all  techniques  by  the  ratio  of  the  average 
received  signal  energy  and  four  times  the  product  of  the  r.  f.  bandwidth,  one-sided 
noise  power  density  and  the  channel  multipath  spread.  The  measurement  times 
are  comparable  for  the  correlation,  multitone  and  the  pulse-pair  techniques.  For 
the  chirp  technique,  the  measurement  time  is  greater  by  a  factor  equal  roughly  to 
the  product  of  the  r.  f.  bandwidth  and  the  multipath  spread. 

The  Correlation  Technique:  The  correlation  technique  for  the  measurement  of 
the  tap-gain  correlation  function  R(S,  f)  or  of  the  scattering  function  Sir^,  f^)  was 
originally  proposed  by  Knilath  (see  Ref.  in  Bello7)  and  analyzed  by  Gallagher 
(see  Ref.  in  Bello  ).  This  technique  involves  the  transmission  of  a  sounding  sig¬ 
nal  with  an  autocorrelation  function  narrow  enough  to  resolve  the  multipath  struc¬ 
tures,  for  example  a  pseudo-random  sequence  of  appropi  iate  length.  At  the 
receiver,  the  received  waveform  is  multiplied  by  a  shifted  replica  of  the  probing 
signal  at  an  offset  frequency  and  the  difference  frequency  component  is  extracted. 
The  autocorrelation  function  of  this  extracted  component  in  the  variable  r  yields 
the  estimate  of  R(£j,  t).  Let  x(t)  be  the  sounding  signal  and  y(t)  be  the  channel  out¬ 
put,  then 

Est(R(?,r)]  =  Cq  J dt  [x(t  -  r)y*<t)|[x*  (t  +  C  -  t)  y(t  +?))  (72) 

and 

y(t)  =  J  x(t  -  t')  g(t,  r  )  dr'  (73) 

where  Cy  is  a  constant. 


59 


7 

The  Multitone  Technique:  The  multitone  technique  is  used  to  estimate  the 
channel  time-frequency  correlation  function  R(D,  t)  by  the  direct  determination  of 
the  cross-correlation  function  between  many  received  tones  spaced  apart  in  fre¬ 
quency.  The  tone  spacing  should  be  larger  than  the  Doppler  spread  and  the  receiver 
filter  bandwidth. 

If  H(ju,  t)  is  the  complex  envelope  of  the  received  tone,  then  an  estimate  of 
the  cross -correlation  function  between  processes  corresponding  to  tones  trans¬ 
mitted  at  frequencies  f  and  f  +  D  Hz  from  the  carrier  frequency  is  given  by 


T/2 

Est[R(n,T)l  =  i  f  H*(ju.t  -  T/2)  H  (ju  +  jn,  t  +  t/2)  dt  .  (74) 

-T/2 


O 

The  Chirp  Technique:  This  technique,  suggested  by  Gallagher  (ref.  in  Bello  ), 
measures  the  channel  time-frequency  correlation  function  R(D,  r).  This  technique 
involves  the  transmission  of  pulses  whose  frequency  varies  linearly  with  time, 
that  is,  a  sequence  of  chirp  pulses  with  different  frequency  slopes.  Each  received 
chirp  pulse  of  a  given  slope  is  multiplied  by  a  replica  generated  at  the  receiver 
but  offset  in  frequency  and  the  difference  frequency  component  is  extracted.  The 
autocorrelation  function  of  this  component  provides  an  estimate  of  R(D,  t)  in  the 
form  of  R((3t,  r)  where  0  is  the  slope  in  Hz  per  second  for  the  chirp  pulse  x(t,  /)). 

Let 


x(t,0)  =  exp  Ijffflt2] 


(75) 


where  Tj  the  pulse  duration,  is  chosen  to  exceed  the  integration  time  T  and  T  must 
be  very  much  larger  than  the  correlation  time  of  the  fading.  It  is  known  that  fad¬ 
ing-correlation  time  is  larger  than  the  path  delay.  Therefore  the  estimate  of 
R((3t,  t)  provided  by  the  chirp-pulse  technique  is  given  by 


T/2 

Est  (R(/)t,  t)l  =  ^  J  x  (t  -  t/2,0)  y*  (t  -  t/2)  x*  (t  +  t/2 ,fi)  y(t  +  t/2)  dt 
-T/2 

(7(5) 

where  y(t)  is  received  complex  envelope  given  by 


y(»)  •=  /  x  (t  -  ?)  g(t,o  d? 


(77) 
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The  Pulse-Pair  Technique:  In  the  pulse-train  probing  technique,  the  tap-gain 
correlation  function  R(?,  t)  is  measured  by  transmitting  a  sequence  (N)  of  pulse 
pairs.  The  separation  t  between  the  members  of  a  pair  is  chosen  larger  than  the 
duration  of  received  pulse  and  the  location  of  the  pulse  pair  is  sufficiently  sepa¬ 
rated  so  that  the  received  signal  consists  of  a  distinguishable  set  of  N  pulse  pairs. 
The  earlier  member  of  the  received  pair  is  delayed  by  t  seconds  and  multiplied  by 

the  conjugate  of  the  later  member.  From  the  resulting  pulse  train,  it  has  been 
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shown  in  Bello  that  the  tap-gain  correlation  function  can  be  obtained. 

2.5  Disciuaion  of  Diversity  and  Adaptivity 

In  this  report  we  consider  two  cases  of  propagation  paths,  one  characterized 
by  a  multipath  spread  L  =  100  microsec,  Doppler  spread  B  =  0. 1  Hz  and  path 
losses  130  dB  (correspondent  to  a  midlatitude  link)  and  another  characterized  by 
multipath  spread  L  =  1500  microsec,  Doppler  spread  B  =  10  Hz  and  path  losses 
150  dB  to  165  dB  (correspondent  to  a  transauroral  link).  We  call  the  first  Case  I 
and  the  second  Case  II.  In  both  cases,  we  have  assumed  that  noise  and  interfer¬ 
ence  reach  overall  levels  of  +30  dB  and  +45  dB  above  KTB. 

The  diversity  situation  is  very  different  in  the  two  cases.  In  Case  I,  the  basic 
diversity  cell  is  a  rectangle  of  size 

[  Xg  =  104  Hz  X  10  sec 

This  means  that  at  a  given  frequency  the  fading  envelope  is  correlated  for  a  time 

interval  of  10  seconds,  while  two  frequencies  simultaneously  emitted  are  correlated 

(fading-wise)  if  their  separation  is  equal  to,  or  less  than  10  KHz,  In  Case  II,  the 

cell  size  is  667  Hz  X  10"1  sec.  To  achieve  diversity,  the  signal's  waveform  must 

occupy  more  than  one  elementary  cell  in  the  diversity  grid.  This  cannot  be 

achieved  by  lengthening  the  duration  of  the  waveform  (it  would  not  be  practical  to 

make  the  length  of  a  bit  several  seconds);  the  only  approach  available  is  to  choose 

the  spectrum  of  the  radiated  waveform  to  be  larger  than  l/L. 

1 8 

By  following  Kennedy  and  Lebow  we  can  see  in  detail  how  to  achieve  the 
required  diversity  in  order  to  assure  an  error  rate  less  than  a  prescribed  maxi¬ 
mum,  for  the  desired  data  rate. 

To  start  with,  we  use  FSK  modulation.  We  will  extend  this  analysis  to  MFSK 
in  a  later  section  of  this  report.  We  alio  ;ssume  that  the  two  transmitted  FSK 
signals  (mark  and  space)  are  composed  of  n  chips,  arranged  either  in  time 


18.  Kennedy,  R.S.  and  Lebow,  I.  L.  (1964)  Signal  desigr  for  dispersive  channels, 
IEEE  Spectrum,  Vol.  1,  pp.  231-237. 


61 


sequence,  in  parallel,  or  in  some  arbitrary  series -parallel  combination.  Let  E 
be  the  total  energy  received  in  all  n  chips.  We  also  assume  that  the  time-band- 
width  product  for  the  radiated  waveform  is  TW  =  1.  The  diversity  per  chip  Z  is 
given  by  the  expressions  contained  in  Table  5.  If  the  mark  (space)  contains  n  chips, 
the  total  signal  diversity  is  Zn. 


Table  5,  Expressions  to  be  Used  in  the 
Computation  of  the  Diversity  Z 


| - ! 

If 

Then 

BT  <  1  and  LW  <  1 

Z  =  1 

BT  <  1  and  LW  >  1 

Z  =  LW 

BT  >  1  and  LW  <  1 

Z  =  BT 

BT  >  1  and  LW  >  1 

Z  =  BL 

In  the  case  of  a  binary  alphabet,  the  probability  of  error  is  computed  by  using 
the  following  expression  (also  plotted  in  Figure  23) 

4#-A(y> 

Pe  s  0.  2  e  *  lN°  (78) 

where  Nq  is  the  noise-power  density  (watts/Hz)  and  the  efficiency  function  A(y)  is 
plotted  in  Figure  24.  The  minimum  for  the  probability  of  error  occurs  when  the 
function  A(y)  is  at  a  maximum.  From  Figure  24,  we  notice  that  a  broad  maximum 
for  A(y)  is  occurring  when  y  a  0.35.  The  corresponding  value  of  A(y)  is 
A(y)  a  0.3. 

By  making  use  of  the  analytical  approach  illustrated  above,  let's  evaluate  how 
much  diversity  we  need  to  counteract  fading  in  the  case  of  a  binary  signalling 
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scheme  meant  to  achieve  a  data  rate  H  =  24  Kilobits/sec  with  an  error  rate  10  . 

We  assume  first  to  be  in  Case  I  path  (L  =  100  microsec,  P  =  0. 1  Hz,  path  losses 
130  dB). 

-4 

From  Figure  25,  in  order  to  achieve  Pe  a  10  ,  we  must  have  E/Nq  A(y)  a  14. 

With  the  best  choice  for  y(ys  0.35)  we  know  that  A(y)  a  0,3,  Therefore,  we 
require 


(0.3)  a  14  ,  #  s  55  (79) 

o  o 


(>2 


9  IP  It  M 


4  *(y) 


Figure  23.  Probability  of  Error 


Figure  24.  Efficiency  Function 
A(y) 


Because,  since 


y  ■  Or,  “  °-35  ■ 


we  conclude  that  we  need 

Zn  =  0.35  X  55  3  20  (81) 

as  total  signal-waveform  diversity.  In  order  to  find  out  how  many  chips  we  need, 
we  must  obtain  from  other  considerations  Z  (diversity  per  chip). 

In  order  to  find  Z  we  must  determine  the  products  BT  and  LAV,  so  that  we  can 
enter  Table  5  and  choose  the  correct  expression  for  Z,  We  have 

R  =  24  Kilobit/sec  and  therefore  T  =  ^  S  41.5  microsec  (82) 

This  is  the  length  of  the  mark  (and  of  the  space)  and  the  bandwidth  that  the  wave¬ 
form.  requires  is 


W  =  ^  =  24  KHz 


We  are  now  in  a  position  to  compute  BT  and  LAV: 


BT  =  0.  1  X  41.  5  X  10  <  1 


G3 


(85) 


LW  =  10“4  X  24  X  103  >  1 


Therefore  from  Table  5,  we  have  that  2  =  LW  and  consequently,  2  =  2.  4.  Because 
then  the  required  total  waveform  diversity  is  Zn  =  20,  we  derive  immediately  how 
many  chips  we  need 

n  =  &  9  (86) 

This  diversity,  required  to  counteract  fading,  can  be  achieved  for  instance  by 
signalling  simultaneously  on  nine  carriers,  each  with  a  24-KHz  bandwidth.  Min¬ 
imum  required  separation  between  two  adjacent  carriers  is  W  +  B  s  24  KHz. 

Counteracting  fading  is  not,  however,  the  only  problem  that  we  have  to  solve 
in  order  to  reliably  transmit  our  24  Kilobit/sec  data  rate.  In  fact,  our  basic  bit 
length  is  41.  5  microsec,  while  the  time  spread  of  the  channel  is  100  microseconds. 
Therefore,  we  cannot  continue  to  signal  in  a  certain  channel,  after  the  transmis¬ 
sion  of  a  bit,  for  at  least  100  microseconds  counted  from  bit's  end,  or  141.  5 
microseconds  from  bit's  beginning.  Consequently,  we  need  a  number  of  groups  of 
frequencies  (each  with  nine  lines)  as  computed  with  the  formula 


T  +  L 
T 


=  1 


100 

41.5 


a  4 


(87) 


We  have  now  all  the  Inputs  required  to  make  an  estimate  of  the  total  bandwidth 
occupancy  required  by  the  link,  as  well  as  the  equipment  complexity: 

Total  bandwidth  occupancy: 

2  X  9  X  4  X 

mark  carriers  groups 

&  per 

space  group 

Total  number  of  carriers  used:  2  X  9  X  4  =  72, 

The  above  figures  apply  to  simplex  communications  and  must  be  multiplied  by 
2  to  achieve  duplex  connection. 

In  order  now  to  complete  the  conceptual  system  design,  we  must  perform  an 
estimate  of  the  required  transmitter  power.  To  do  this,  we  compute  first  the 
expected  noise  and  interference  level  in  reception.  As  mentioned  in  previous 
sections,  we  have  adopted  at  HF  two  values  for  this  level:  +30  dB  and  +45  dB  above 
KTB.  Therefore,  for  the  first  alternative,  we  have 


24  KHz  =  1.728  MHz 

bandwidth 

per 

carrier 
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(88) 


Nq  =  103(KT)  =  103  1.38  10'23  3  102  =  4.  14  10"18  watts /Hz  . 

The  required  signal  energy  is  therefore 

E  =  55  NQ  =  2.  28  10"16  joules  (89) 

Consequently,  the  required  signal  power  is 

Pr  =  jjr  =5.48  10"12  watts  .  (90) 

Because  the  path  losses  for  Case  I  are  130  dB,  inclusive  of  antenna  gains,  the 
transmitter  must  have  a  power  level 

PT  =  1013  Pr  =  55  watts  .  (91) 

Under  the  assumption  that  noise  and  Interference  amount  to  +45  dB  above  KTB,  the 
required  level  of  transmitter  power  is 


PT  =1.74  K  watt  .  (92) 

We  have  to  see  now  how  the  situation  changes  when  the  propagation  path  becomes 
the  one  of  Case  II  (L  =  1500  microsec,  B  =  10  Bz,  path  losses  150  dB  to  165  dB). 
Those  are  typical  conditions  foi  transauroral  propagation.  The  diversity  conditions 
change  radically.  We  now  have 

BT  =  10  41.5  10-G  <  1 

(93) 

LW  =  1.  5  10'3  24  103  .»  1 


and  the  chip  diversity  is  now  Z  =  LW  =  24  1.  5  =  36.  Table  6  illustrates  the  trans¬ 
mitter  power  requirements,  in  this  case. 

Because  the  total  diversity  requirement  for  the  waveform  is  Z  n  =  20,  a  single 
chip  (n  =  1)  is  more  than  enough.  The  number  of  carriers  per  group  is  decreased, 
when  going  from  Case  I  to  Case  II,  from  9  to  1;  however,  there  is  a  severe  deteri¬ 
oration  in  the  situation  of  inter-symbol  interference.  In  fact,  we  need  now  a  much 
larger  number  of  groups: 
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Table  6.  Transmitter  Power  Requirements  in  Transauroral 
Paths  for  Binary  Waveforms 


1 

Noise  and  Interference 

Required  Transmitter 

Path  Losses 

Level 

Power 

150  dB 

+30  dB  above  KTB 

5,  5  K  watts  1 

+45  dB  above  KTB 

I'M  K  watts  j 

1G5  dB 

+30  dB  above  KTB 

174  K  watts 

+45  dB  above  KTB 

5,  5  M  watts 

T  +  L 
T 


=  1 


+  1500  „ 

+  TTF  238 


(94) 


This  is  because  we  must  wait  1500  mlcrosee  from  the  end  of  a  certain  bit  before 
being  able  to  go  back  to  signal  in  the  same  channel. 

The  estimate  of  total  bandwidth  occupancy  and  of  equipment  complexity  now 
yields: 

Total  bandwidth  occupancy; 

2X1  X  38 

mark  carriers  groups 

&  per 

space  group 

Total  number  of  carriers  used:  2  X  1  X  38  -  70 

Not  withstanding  the  large  differences  between  Case  1  and  Case  11,  as  far  as 
the  diversity  condition  is  concerned,  there  is  a  close  similarity  between  the  two 
cases  in  terms  of  total  bandwidth  occupancy  and  total  number  of  carriers  used. 
Actually,  this  is  not  totally  unexpected  because  going  from  a  path  with  little  time 
spread  to  one  with  large  multipath  spread,  we  gain  in  terms  of  protection  against 
fading  (for  a  given  waveform),  but  we  lose  m  terms  of  inter-symbol  interference. 
These  results  give  a  clear  indication  how  the  adaptivity  of  the  system  could 
work.  If  we  have  available  for  instance  70  spectral  lines  in  the  IIF  band,  each  with 
a  bandwidth  of  24  KHz,  we  could  arrange  72  of  them  into  four  groups  of  nine  car¬ 
riers  when  the  path  is  Case  I  type  or  \\°  could  use  all  7G  of  them  (11)  u,  38  groups, 
each  with  a  single  cat  rier.  This  adaptive  adjustment  to  path  conditions  could  be 
made  based  on  the  results  of  path  sounding,  channel  probing,  as  well  as  oil  the 


X  24  KHz  =  1.82  MHz 

bandwidth 

per 

carrier 


GG 


basis  of  noise  and  interference  measurements  performed  for  each  frequency  in 
use.  A  control  unit  would  use  this  information  to  decide  on  the  best  possible 
arrangements  of  frequencies  and  groups,  for  local  use  at  the  terminal  where  these 
determinations  are  made  and  for  transmission  to  the  other  terminal,  to  adapt  the 
parameters  of  both  station  equipments  to  the  conditions  of  the  path. 

2.6  The  Use  of  M-ary  Codes 

A  substantial  improvement  in  the  link's  performance  can  be  obtained  by  using 
M-ary  codes.  When  the  FSK  modulation  is  used,  the  corresponding  approach  is 
called  MFSK.  With  these  codes,  it  is  possible  in  principle  to  transmit  more 
information  through  the  channel  per  unit  time.  The  M-ary  waveform  considered 
by  us  consists  in  the  transmission  of  one  out  of  M  frequencies,  equally  likely,  so 
that  the  number  of  bits  thus  transmitted  is  logg  M.  It  would  be  logg  2  =  1  if  the 
transmission  were  by  binary  alphabet. 

In  our  study,  we  have  considered  a  code  with  M  =  8.  The  adoption  of  this  code 
makes  it  possible  to  achieve  a  data  rate  of  if-'  Kiloblt/sec  with  a  pulse  length 
T  =  41,  5  X  log0  8  =  125  microsec,  corresponding  to  a  bandwidth  W  -  8  KHz.  With 
M-ary  codes,  the  probability  of  error  is  worse  than  for  the  binary  case.  In  fact 

pbit,  M  =  ok  -  1  (M  "  l>  Pbit-  Binai7  (95) 


where  k  =  log,  M. 

*• 

If  we  want  to  keep  ^  =  10  as  in  the  binary  case,  we  must  have,  for 
M  =  8 


P 


bit,  Binary 


2.5  10'5 


(98) 


From  Figure  23,  we  have  that  this  probability  of  error  (binary-equivalent)  is 
obtainable  with  E/Nq  A(y)  =  17.  The  optimum  binary-equivalent  E/Nq  ratio  occurs 
for  A(y)  =0.3  and  can  be  derived  from  the  usual  equation 


We  obtain 


(97) 


67 


The  diversity  Z  per  chip  can  be  computed  as  follows 
Case  I 

LW  =  0. 8  <  1 

BT  =  125  10"7  <  1 

therefore,  Z  =  1 
Case  II 

LW  *  1.5  10"3  8  102  =  12  >  1 
BT  *  125  10*G  10  =  12.5  10"4  <  1 


(98) 


(99) 


therefore,  Z  =  LW  s  12. 

We  conclude  that  in  Case  I,  because  the  required  total  diversity  ig  Z  n  =  20, 
we  need  only  20  carriers,  each  with  a  bandwidth  of  8  KHz.  In  Case  II,  a  total  of 
two  carriers  is  more  than  enough  to  counteract  path  fading.  In  order  to  counteract 
inter-symbol  interference,  the  number  of  groups  that  we  require  can  be  estimated 
as  follows 

Case  I 


(100) 


Case  II 


T  +  L 
T 


a  13 


(101) 


We  have  now  all  the  elements  to  make  an  estimate  of  the  total  bandwidth  occupancy 
required  by  the  MFSK  (with  M  =  8)  link,  as  well  as  of  the  equipment  complexity 


Case  I  —  Total  bandwidth  occupancy  (for  simplex) 

2  X  8  X  20  X  2  X  8  KHz  =  5.  12  MHz 

mark  M-ary  carriers  groups  bandwidth 
&  Code  per  per 

space  group  carrier 

Total  number  of  carriers  used:  2X8X  20  X2  =  640  (for  simplex) 
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or  (104) 

E  =  2.  5  10"15  joules  . 

The  required  levels  of  transmitter  power  are  given  in  Table  7.  It  can  be  seen 
that  there  is  a  factor  of  almost  10  improvement  with  respect  to  the  binary-system 
approach,  at  the  expense  of  bandwidth  occupancy  and  equipment  complexity. 

The  power  densities  (watts/Hz)  radiated  by  the  proposed  M-ary  coded  emis¬ 
sions  are  given  in  the  fourth  column  of  Table  7.  The  equivalent  power  of  a  trans¬ 
mitter  that  wouiu  generate  equal  power  density  in  a  3-KHz  voice  channel  is  given 
in  the  fifth  column  of  the  same  table.  These  levels  are  relatively  low,  with  the 
exception  of  the  one  on  the  last  line.  We  have  also  to  point  out  that  our  transmitter 
would  be  even  less  bothersome  because  the  emissions  would  not  stay  consistently 
on  the  same  frequencies,  but  would  wander  around  to  follow  the  adaptivity  instruc¬ 
tions  generated  by  the  microprocessors  and  control  logics. 
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Table  7.  Transmitter  Power  Requirements  in  HF  Paths  for  M-ary 
Transmission  With  M  =  8 


Noise  and 

Required 

Radiated 

Equivalent 

Interference 

Transmitter 

Power 

Power  in  a 

Path  Losses 

Level 

Power 

Density 

3 -KIIz  Channel 

Midlatitude 

Path 

+30  dB  above  KTB 

6.  67  watts 

1.3  10'°  watts 

4  milliwatts 

130  dB 

+45  dB  above  KTB 

211  watts 

4.  11  10‘5"H* 

123  milliwatts 

Transauroral 

Path 

150  dB 

+30  dB  above  KTB 

+45  dB  above  KTB 

667  watts 

21,1  K  watts 

2.4  10"4  watts 
6.34  10"3  Hz 

0.  6  watts 

19  watts 

105  dB 

+30  dB  above  KTB 

21.  1  K  watts 

6.34  10'3 

2  10"1 

19  watts 

+45  dB  abvve  KTB 

667  K  watts 

600  watts 

Note:  Pulsowidth  •  125  microsec 


2.7  Nan-Adaptive  One-Way  Link  for  Information  Transfer  Without 
Feed-Hack  Acknowledgement 

There  are  eases  in  which  it  is  neither  operationally  possible  nor  advisnble  to 
establish  a  two-way  link  with  feed-back  acknowledgement  provisions  between  two 
terminals.  In  such  cases,  adaptivity  cannot  be  embodied  into  the  system  nnd  a 
substantial  decrease  of  communications  reliability  may  result.  This  is  the  case 
that  was  called  Type  B  Link  in  Section  2. 1. 

The  way  that  this  case  is  handled  traditionally  is  to  repeat  the  message  on 
several  frequencies.  We  present  in  this  section  an  alternative  method  based  on 
the  use  of  M-ary  transmission  with  a  large  value  of  M.  Because  we  have  already 
introduced  the  M-ary  transmission  method  in  Section  2,  6  and  in  Appendix  A,  we 
limit  ourselves  here  to  a  numerical  example. 

We  assume  io  have  M  =  G4  and  we  still  keep  the  desired  data  rate  at  24  Kilobit/ 
sec  and  the  specified  error  rate  at  10~4.  Pulse  length  is  now  T  =  41.5.  log2  64  = 
249  microsec,  that  corresponds  to  a  bandwidth  of  4  KIIz. 

We  use  the  same  equation  for  the  probability  of  error  that  was  used  in  Section 

2.  6 


gk-i 

^bit,  M  9k  _  j  1  bit,  binary 


(105) 
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where,  as  usual,  k  =  logg  M.  In  our  case  we  have  therefore  that  k  =  6,  and 

pbtt,  binary  s  TT  "  3  W’6  <‘°6> 

The  function  binary  is  re^ated  to  E/N0  ratio  ^  formula 

pblt,  binary  1  Pe  =  °'2  “I1  K  <E/N0>  A(y>}  (107) 

where  the  symbols  have  the  same  meaning  as  indicated  in  Section  2.  5. 

From  the  formula  above,  we  obtain,  by  adopting  as  usual  A(y)  =  0.  3 

E/N0  =  73.77  (108) 

Therefore,  Z  n  =  0.  35  E/Nq  *  26.  For  the  midlatitude  path  (Case  I),  Z  =  1  (be¬ 
cause  LW  <  1  and  BT  <  1),  hence  n  =  26.  For  the  transauroral  path  (Case  II,  we 
have  Z  *  LW  =  6  (because  nov  LW  >  1,  BT  <  1),  hence  n  =  5. 

In  order  to  counteract  the  inter-symbol  interference  problem,  we  have  to 
compute 

Case  I 

=  2  (109) 

Case  II 

»  7  (110) 

With  this,  we  now  have  all  the  parameters  necessary  for  the  evaluation  of  the 
required  bandwidth  occupancy  and  of  the  number  necessary  for  the  carriers: 

Case  I  —  Total  bandwidth  occupancy  (f<  r  simplex) 

2X  64  X  26  X2X4  KHz  =  26.  62  MHz 
Total  number  of  carriers  used  (for  simplex) 

2  X  64  X  20  X  2  «  6656  carriers 

As  the  reader  can  see,  this  numerical  example  was  worked  out  to  correspond 
to  the  total  bandwidth  occupancy  of  the  entire  HF  band  considered  in  our  study: 
from  3  MHz  to  30  MHz  (27-MIIz  band  considered  as  available  in  p.  inciple).  This 
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was  a  consequence  of  the  adoption  of  the  value  M  =  64.  A  calculation  that  would  be 
of  Interest  to  perform  in  follow-up  activity  would  be  the  estimate  of  the  increase 
in  error  rate  when  portions  of  the  HF  spectrum  (3  to  30  MHz)  disappear  in  propaga¬ 
tion  outages. 

For  the  transauroral  paths,  we  have 
Case  II  —  Total  bandwidth  occupancy  {for  simplex) 

2X64X5X7X4  KHz  =  17.  92  MHz 
Total  number  of  carriers  used  (for  simplex) 

2  X  64  X  5  X  7  =  4480  carriers 


As  far  as  chip's  energy  is  concerned,  we  have 


E/N0  •  *  12.3  (111) 

in  the  binary-equivalent  system.  Because  the  values  that  we  have  adopted  for  Nq 
are  as  follows: 


Nq  =  4. 14  10" 18  watts/Hz  (+30  dB  above  KT) 
or 

Nq  =  1.3k  i0'lR  watts/Hz  (+45  dB  above  KT) 
we  obtain  the  required  signal  energy  levels 
12  =  50.9  10'1!i  joules 


(112) 


or 


(113) 


E  =  16.  11  I0"l(i  joules  . 


Table  8  provides  the  values  of  the  transmitter  power  levels  necessary  to 
achieve  the  signal  energies  above  (oulse  width  =  249  microsoc). 
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Table  8.  Transmitter  Power  Requirements  in  HF  Paths  for  Non-Adaptive 
M-ary  Transmission  With  M  =  64 


Radiated 

Noise  and 

Required 

Power 

Equivalent 

Interference 

Transmitter 

Density 

Power  in  a 

Path  Losses 

Level 

Power 

(watts/Iiz) 

3 -KHz  Channel 

Midlatitude 

Paths 

0.075  10"(> 
2.419  10"G 

130  dB 

^30  dB  above  KTB 

2  watts 

0.  225  milliwatts 

♦■45  dB  above  KTB 

(54,  4  watts 

7. 258  milliwatts 

Transauroral 

Paths 

♦30  dB  above  KTB 

200  watts 

11.  1G  I0'(i 

33.  5  milliwatts 

150  dll 

+45  dB  above  KTB 

C,  44  kilowatt 

0.359  10‘3 

1. 1  watts 

1(55  dB 

♦  30  dB  above  KTB 

(5. 44  kilowatt 

0.359  10~3 
11.33  10'3 

1.  1  watts 

♦45  dB  above  KTB 

203.  (!  kilowatt 

34  watts 

Note:  Pulsewidth  249  microsec 


3.  I'llOI’OSKI)  KXHKRIMBNTAI.  ACTIVITY 

Experiments  on  the  diversity  and  adaptivity  aspects  of  the  proposed  scheme 
could  be  performed  by  adopting  an  evolutionary  approach  that  could  start  from  the 
analysis  of  existing  digisonde  data  (records  of  the  IF  output  are  essential  in  order 
to  perform  measurements  of  the  Doppler  spread),  and  from  a  review  of  the  state- 
of-the-art  of  the  required  instrumentation  (frequency-agile  transmitters  and 
receivers,  broadband  HF  antennas  and  matching  units,  high-speed  frequency 
switching,  etc.). 

Later  on  we  could  construct  the  breadboard  of  a  one-way  link  designed  accord¬ 
ing  to  the  principles  of  the  system  approach  illustrated  herewith  and  we  could  ob¬ 
tain  valuable  data  of  the  "path -sounding"  and  "channel-probing"  category.  These 
data  would  provide  a  reliable  picture  of  the  diversity  aspects  of  midlatitude  and 
transauroral  paths. 

Fundamental  parameters  that  must  be  gathered  both  for  the  midlatitude  and  the 
transauroral  case  are  the  following: 

1.  Statistics  of  path  losses,  to  be  summarized  in  a  diagram  "percentage  of 
time  that  the  path  losses  are  lower  than  the  abscissa"; 

2.  Statistics  of  multipath  spread,  again  leading  to  the  diagram  "percentage  of 
time  that  the  multipath  spread  is  smaller  than  the  abscissa"; 
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3.  Statistics  of  Doppler  spread,  again  to  be  summarized  in  the  diagram 
"percentage  of  time  that  the  Doppler  spread  is  smaller  than  the  abscissa"; 

4.  Statistics  of  the  noise  and  interference  levels,  leading  to  the  diagram: 
"percentage  of  time  that  the  level  of  noise  plus  interference  is  smaller  than  the 
abscissa. 

The  measurements  above  must  be  carried  out  for  extended  periods  both  for  the 
midlatitude  and  the  transauroral  case,  and  must  be  performed  at  a  number  of  fre¬ 
quencies  in  the  band  3  MHz  to  30  MHz  at  least  equal  to  the  number  of  frequencies 
mentioned  in  Table  4.  In  principle,  the  separation  between  two  adjacent  frequencies 
should  actually  be  smaller  than  or  equal  to  the  coherent  bandwidth  of  the  path  (in 
the  models  adopted  in  Section  2,  this  is  100  KHz  for  the  midlatitude  path  and  only 
666  Hz  for  the  transauroral  path).  However,  the  adoption  of  this  criterion  would 
lead  to  rather  cumbersome  instrumentation,  and  our  suggestion  is  to  start  first 
with  a  number  of  frequencies  as  indicated  in  the  Table  4  mentioned  above.  This 
would  be  sufficient  for  the  limited  albeit  Important  scope  of  exploring  the  feasibil¬ 
ity  of  the  specific  scheme  proposed  in  Section  2,  where  the  closest  that  the  spot 
frequencies  were  allocated  was  8  KHz. 

Once  that  the  one-way  link  tests  would  be  completed,  we  could  perform  exper¬ 
iments  with  an  expanded  breadboard,  that  could  make  it  possible  to  test  the  adaptiv¬ 
ity  features  of  the  proposed  scheme.  This  would,  of  course,  require  the  two-way 
link  between  the  two  terminals. 


4.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  conclusions  of  our  study  on  adaptive  utilization  of  11F  propagation  paths 
Indicate  that  by  spreading  the  radiated  power  across  several  MHz  of  bandwidth 
occupancy,  it  is  possible  to  counteract  effectively  both  waveform  fading  and 
inter-symbol  interference.  By  this  method,  data  rates  of  24  Kbit/sec  with  10"4 
error  rates  are  possible.  If  this  is  confirmed  by  the  proposed  experimentation, 
HF  would  acquire  a  reliability  and  channel  capacity  of  the  degree  that  is  enjoyed 
by  other  communications  media.  The  penalty  that  must  be  paid  to  achieve  these 
results  is  equipment  complexity  and  bandwidth  occupancy.  Concerning  the  first 
point,  modern  advances  in  microprocessor  technology,  high-density  packaging, 
frequency  agility,  etc.,  offer  concrete  promises.  As  far  as  the  second  point  is 
concerned,  the  power  of  the  transmitter  is  so  spread  that  the  link  would  hardly 
interfere  (and  for  not  very  long  time  intervals,  because  of  frequency  wandering) 
with  a  receiving  site. 
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Our  recommendation  is  that  the  present  study  be  followed  by  additional  R&D 
activity  in  three  basic  directions:  a  deeper  understanding  of  the  propagation 
properties  of  the  two  types  of  paths  considered,  a  gathering  of  experimental  data, 
and  an  engineering  study  of  the  availability  and  applicability  of  such  modern  tech¬ 
nological  breakthroughs  as  microprocessors,  frequency-agile  transmitters  and 
receivers,  broadband  HF  antennas  and  matching  units,  switch  circuits,  decision 
logics,  etc. 
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Appendix  A 

M  ary  Transmiislon  Methods 

Al.  GENERA!. 

In  this  appendix,  the  p<.  rformance  of  three  basic  data-transmlsslon  methods, 
namely,  m-level  Phase  Shi:!  Keyed,  m-level  Amplitude  Shift  Keyed  and  m-level 
(orthogonal)  Frequency  Shill  Keyed  are  compared.  The  problem  addressed  is  one 
of  detecting  a  known  signal  in  a  white  noise  background  with  minimum  probability 
of  error.  The  approach  adopted  is  that  of  Arthurs  and  Dym^1  who  address  this 
problem  from  a  geometric  point  of  view. 

The  analysis  of  data -transmission  systems  is  commonly  based  on  the  following 
model.  There  is  assumed  to  exist  a  message  source  generating  a  stream  of  equally 
likely  messages,  Mj,  Mg,  ....  M|n,  into  a  waveform  generator  having  available 
an  alphabet  of  m  distinct  waveforms,  Sj(t),  Sg(t>,  ...,  Sm(t),  each  of  duration  T 
(and  necessarily  finite  energy).  On r  waveform  is  transmitted  every  T  seconds,  the 
choice  of  waveform  depending  in  some  fashion  "i  the  incoming  message  and  possi¬ 
bly  on  the  waveforms  transmitted  in  preceding  time  slots.  The  medium  coupling 
the  transmitter  to  the  receiver  is  assumed  to  add  stationary  white-zero  mean- 
Gaussian  noise  to  the  transmitted  signal  but  otherwise  is  assumed  to  be  distortion 
free.  It  is  generally  further  assumed  that  the  receiver  is  time-synchronized  with 


Al.  Arthurs,  E.  and  Dym,  11.  (19G2)  On  the  optimum  detection  of  digital  signals 
m  the  presence  of  white  Gaussian  noise,  IRE  Trans,  on  Communications 
Systems,  Vol.  CS-10,  pp.  336-372. 
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the  transmitter  (synchronous  detection).  Sometimes  it  is  also  assumed  that  the 
receiver  is  phase-locked  to  the  transmitter  (coherent  detection).  Here  we  shall 
always  assume  time  synchronism  but  shall  distinguish  between  coherent  and  inco¬ 
herent  detection. 

The  problem  we  are  generally  interested  in  solving,  given  this  model  (see 
Figure  Al),  is  how  to  design  the  receiver  so  that  it  makes  as  few  errors  as  possi¬ 
ble.  Furthermore,  assuming  than  an  optimum  receiver  (optimum  in  the  sense 
that  it  will  make  fewer  errors  in  the  long  run  than  any  other  receiver)  is  construc¬ 
ted,  we  are  interested  in  calculating  its  error  rate, 

TRANSMITTER  CHANNEL  RECEIVER 


In  this  model,  following  assumptions  are  made: 

1.  Each  signal  waveform  is  transmitted  with  equal  probability. 

2.  The  transmitter  is  subject  to  an  average  power  limitation,  E/T  (watts), 
where  T  is  the  duration  of  each  transmitted  signal  waveform. 

3.  The  received  signal  is  the  sum  of  the  transmitted  signal  and  stationary 
white-zero  mean-Gaussian  with  double-sided  spectral  density  NQ  (watts/Ilz). 

4.  The  receiver  is  in  time  synchronism  with  the  transmitter;  however,  the 
distinction  is  made  whether  the  receiver  is  phase-locked  to  the  transmitter. 

5.  The  received  signal  is  processed  by  a  maximum-likelihood  detector  except 
in  the  ASK  incoherent  case. 

6.  In  FSK  and  PSK  cases,  the  transmitted  sinusoidal  pulses  contain  equal 
energy  E  whereas  in  the  ASK  case  the  amplitudes  (square  root  of  the  energy)  cf 
the  transmitted  pulses  are  uniformly  spaced  starting  with  zero.  Furthermore, 

FSK  signals  are  orthogonal. 
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A2.  DISCUSSION 


A2.1  Introductory  Remarks 

In  the  following  discussions,  we  will  compare  the  performance  of  the  M-ary 
data -transmission  systems,  when  a  fixed-waveform  duration  is  transmitted,  and 
when  waveform  duration  is  varied  to  maintain  a  fixed-signalling  rate.  In  Sections 
A2.2  and  A2.3  we  will  provide  the  curves  (only  for  FSK  case)  of  the  probability  of 
error  (10  log jq  pg)  as  a  function  of  signal-to-noise  ratio  (in  decibels),  with  the 
number  of  levels  m  appearing  as  a  parameter.  In  Section  A2.4  we  compare  these 
systems  on  the  basis  of  bandwidth. 

A2.2  Companion*  Aisuming  Fixed-Waveform  Duration 

It  is  known  that  increasing  m  (the  number  of  waveforms  transmitted  in  anytime 
T)  tends  to  increase  the  probability  of  error  whereas  increasing  the  energy  content 
in  each  transmitted  signal  (that  is,  the  signal-to-neise  ratio)  tends  to  decrease  the 
probability  of  error.  Increasing  m  introduces  the  most  degradation  in  the  ASK  ease, 
somewhat  less  degradation  in  the  PSK  case  and  comparatively  little  degradation  m 
the  FSK  case  (Figures  A2  and  A3).  In  fact,  for  large  m,  the  average  probability  of 
error  of  an  FSK  system  should  be  smaller  than  that  of  a  PSK  or  an  ASK  system, 

If  one  compares  these  systems  on  the  basis  of  coherency  vs.  incoherency,  the  PSK 
system  suffers  the  most  degradation  in  performance  due  to  lack  of  coherence. 
Calculations  show  that  for  large  m,  the  cost  of  incoherence  is  a  3-dU  degradation 
in  signal-to-noise  ratio.  For  small  m,  ‘.he  degradation  is  somewhat  less.  In  the 

ASK  and  FSK  cases  (Figure  A4)  it  has  been  shown  that  for  a  small  probability  of 

•5 

error  (say,  pe  <  10  )  the  degradation  between  coherent  and  incoherent  is  of  the 

order  of  a  decibel  or  less. 

A2.3  Cumpariion*  Aiuitmmg  a  Fixed-Signaling  Kate 

Instead  of  maintaining  the  duration  of  each  transmitted  pulse  at  a  fixed  value 
T  independent  of  nt,  we  now  consider  the  equally  valid  constraint  of  a  fixed  rate  R 
(we  are  still  assuming  that  the  transmitter  is  aver  age -power  limited).  Under  tins 
constraint,  we  can  allow  a  longer  time  duration  for  each  waveform  tn  a  multilevel 
system  and,  hence,  increase  the  energy  content  of  the  transmitted  signals.  Assum¬ 
ing  that  the  rate  at  which  data  is  being  transmitted, 

log„  m 

R  =  — -  (bits/sec)  ,  (At) 

m 
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Figure  A2.  Probability  of  Waveform 
Error  (r  .-Level  FSK  Coherent)  Assum¬ 
ing  That  the  Duration  of  Each  Signal  Is 
Fixed  Independently  of  m 


Figure  A3.  Probability  of  Waveform 
Error  (m-Level  FSK  Incoherent)  As¬ 
suming  That  the  Duration  of  Each  Signal 
Is  Fixed  Independently  of  m 


is  maintained  constant,  it  follows  that 


(A  2) 


and,  therefore,  that 

Tm  =  T2  log2  m  (A3) 

since  the  transmitter  is  average-power  limited. 


E. 

T 


m 


m 


(A4) 


Combining  Eqs.  (A3)  and  (A4)  yields 


E 


m 


Eg  log9  m 


(A  5) 
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FigureA4.  Comparison  of  Col  .--rent  Figure  A5.  Probability  of  Waveform 
FSK  and  Incoherent  FSK  Error  (m-Level  FSK  Incoherent)  As¬ 

suming  That  Signal  Duration  Is  Adjusted 
to  Keep  the  Data  Hate  Constant 


from  which  it  follows  that 


10  logl0 


E2 

10  log  TjT^-  +  10  lotf10  (log2  m) 
“  o 


(AG) 


That  is  to  say,  under  the  assumptions  of  constant  rate  an  m-level  system  has  a 
signal-to-ncise  ratio  advantage  of  10  log1Q  (log2  m)  dB  over  its  two-level  counter¬ 
part,  If  one  computes  me  probability  of  error  for  these  cases,  one  finds  that  the 
multilevel  PSK  and  ASK  systems  are  inferior  in  performance  to  their  two-level 
counterparts  whereas  FSK  systems  seem  to  improve  in  performance  with  increas¬ 
ing  m  (see  Figure  A 5). 


A2.4  Discussion  of  Bandwidth 

Thus  far,  all  comparisons  have  been  made  under  the  assumption  that  a  distor¬ 
tionless  "wideband"  Gaussian  channel  is  available.  In  practice,  ihe  bandwidth 
allotted  to  any  one  transmitter  is  generally  limited  and,  hence,  the  relative  effi¬ 
ciency  with  which  it  uses  the  available  band  .  idth  of  prime  interest.  The  parameter 
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r  is  defined  as  the  ratio  of  the  rate  at  which  information  is  being  transmitted, 

R  =  log2  m/T  bits/sec,  to  the  Nyquist  ■  *te  of  transmission,  2B  bits/sec,  that  is 


r  = 


R  l°S2  m 


2B 


2BT 


(A  7) 


It  is  quite  difficult  to  define  B  precisely;  for  the  purposes  of  the  present  discussion, 
if  a  sinusoid  of  duration  T  and  frequency  f  will  be  passed  with  negligible  distortion 
by  an  ideal  filter  with  passband  1.  5/T  centered  around  fQ,  It  follows  therefore, 
that  for  the  PSK  and  ASK  modulation  schemes  wherein  the  frequency  of  the  pulses 
sent  in  each  time  slot  is  fixed  the  required  transmitter  bandwidth  B  a  1,  5/T  and, 
correspondingly,  the  bandwidth  efficiency 


r  = 


log2 


m 


3 


(A  8) 


in  the  FSK  case,  however,  assuming  a  separation  of  1/T  between  adjacent 
tones,  an  m-level  transmitter  requires  a  bandwidth 


B 


>n  +  0. 5 
T 


(A  9) 


and  which  case 
log2  m 


r  * 


2m  +  1 


(A  10) 


In  the  FSK  coherent  case,  however,  adjacent  signals  need  only  be  separated 
by  a  frequency  difference  of  1/2T  to  maintain  orthogonality  and,  hence,  the  required 
bandwidth  may  be  reduced  to 


B  = 


m  +  2 
2T 


(All) 


resulting  ha  a  bandwidth  efficiency  of 


r  = 


log2  m 
^Tm 


(A  12) 
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It  is  apparent  from  Figure  A 6  that 
simple  multilevel  orthogonal  FSK  sys¬ 
tems,  even  under  idealized  operating 
conditions,  are  inefficient  users  of 
bandwidth.  Physically,  the  reason  is 
clear;  in  any  time  period  T  only  a 
fraction  of  the  total  system  bandwidth, 
namely,  that  occupied  by  the  particular 
tone  transmitted,  is  utilized.  Thus, 
we  see  that  although  the  number  of 
levels  of  an  FSK  system  may  be  in¬ 
creased  with  relatively  little  degrada¬ 
tion  in  performance,  there  is  corre 
spondingly  an  Increase  in  the  bandwidth 
required  by  the  system  to  operate.  In 
contrast,  multilevel  ASK  and  PSK  mod¬ 
ulation  schemes  utilize  the  bandwidth 
more  efficiently  but  the  probability  of 
error  increases  with  m. 

A2.5  Signal  Selection  for  M  ary  Trarwmi«»ion 

For  the  complete  evaluation  of  a  communi  ation  system  and  the  related  coding 
schemes,  one  needs  to  determine  not  only  the  error  probability  as  a  function  of  the 
signal-to-noise  ratio  and  data  rate,  but  also  the  bandwidth  occupancy  of  the  trans¬ 
mission  channel.  The  most  direct  definition  of  bandwidth  occupancy  of  a  given 
channel  is  the  minimum  frequency  separation  between  channels  such  that  signals 
of  one  channel  have  no  effect  on  the  coder  of  any  other.  Let  a  signal  of  the  «-th 
channel  be  denoted  by  x^^(t)  and  one  of  another  channel  x^^(t).  The  o-th  channel 
decoder  operates  on  any  received  signal  y(t)  with  the  M  possible  transmitted  signals 
x^^(t)  of  the  given  channel  to  form  the  quantities 

'I' 

T'  /  y(i)  xft(l)(t)  dt  .  (A  13 1 

o 

Thus,  there  will  be  no  effect  on  the  decoder  of  either  channel  due  to  ihe  signal  of 
the  adjacent  channel  provided 


NUMBER  OF  LEVELS 

Figure  AG.  Plot  of  Bandwidth  Efficiency 
ns  a  Function  of  the  Number  of  Levels 
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f  x  (l)(t)  xJi\t)  dt  =  0  (A  14) 

J  a  P 
o 

th 

for  all  i  and  j  and  for  all  a  and  p.  Denoting  the  elementary  signals  of  the  a  chan¬ 
nel  by  SQ  and  S^,  and  those  of  the  0^  by  and  S^,  the  condition  becomes 

T/n 

/  Sr(t)  Sm(t)  dt  =  0  (A15) 

0 

for  r  and  m  +  0,  1,  where  each  output  signal  x^(t)  of  duration  T  seconds  consists 
of  a  sequence  of  n  signals,  each  of  duration  T/n  seconds,  chosen  randomly  from 
the  two  elementary  signals  SQ(t)  and  S^t).  Thus  the  bandwidth  occupied  by  x^(t) 
depends  on  the  time  functions  SQ(t)  and  Sj(t), 

The  above  orthogonal  condition  is  satisfied  if  the  elementary  signals  of  the 
o-th  channel  are  taken  to  be 

Sr(t)  -  •/IS'cos  (^?)t  (A16) 

those  of  the  other  channel  are 


s|.(t)  =  -  v/lTS' cos  (£j2)  t 


(A  17) 


Thus  it  is  clear  that  the  minimum  frequency  separation  or  the  bandwidth  occupancy 
of  each  channel  is  (ffn/T)(rad/sec]  or  (n/2T)  |  Hz].  Other  channels  can  obviously 
be  placed  at  radian  frequencies  which  are  multiplies  of  »n/T  away. 

The  transmitted  signals  are  readily  generated  by  coding  each  elementary  sig¬ 
nal  which  itself  is  generated  cither  by  balanced  modulating  the  carrier  or  by  phase 

modulating  it  by  ±90°.  At  this  point  it  is  useful  to  review’  and  discuss  briefly  the 

K  (i) 

basic  physical  parameters.  A  set  of  M  =  2  equiprobable  signals  x'  ;(t),  each  of 
duration  T  seconds  are  received  nt  a  power  level  of  S  watts  with  additive  stationary 
white  Gaussian  noise  of  spectral  density  N  watts/Hz.  The  relationships  among 
signals  which  are  pertinent  to  the  reception  problem  are  their  normalized  inner 
products 


J  x(0(t)  x(^(t)  dt 

Q 

py  '  "  st 


(a  is: 
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The  error  probability  for  the  optimal  decoder  is  a  function  of  ST/N  and  the 
[p.j]  matrix.  When  x'  '(t)  is  generated  as  a  sequence  of  n  elementary  binary  sig¬ 
nals  p . .  is  expressible  in  terms  of  the  number  of  agreements  and  disagreements 
lJ  (!)  (l\ 

in  the  binary-code  vectors  from  which  x  (t)  and  x  J  (t)  was  generated.  For,  in 

this  case 

n  r(T/n) 

Ptj  =  ^  Z  /  x(i)(t>  x(j)(t)  dt  .  (Ali 

r=l  (r-l)T/n 

For  the  elementary  signals,  each  term  of  the  summation  is  either  +ST/n  or  -ST/n 
depending  on  whether  the  code  symbols  for  and  x^  during  the  specific  sub- 
interval  were  the  same  or  different.  Hence  for  binary  signals, 

#  agreements  -  #  disagreements 


Any  one  of  the  M  equiprobable  transmitted  signals  contains  logg  M  bits  of  informa 
tlon.  Since  each  signal  requires  T  seconds  for  transmission,  the  transmitted 
information  rate  is  given  by 


log9  M  K 

R  s  =  £  bits/sec 


The  time  required  to  transmit  one  bit  is 


/3  log2M 


7m  =  l  sec/bit  • 


In  order  to  compare  communication  systems  properly  for  different  values,  the 
error  probability  should  be  considered  as  a  function  of  ST^/N'o,  which  is  the 
received  signal  energy  per  bit,  divided  by  the  noise  spectral  density.  This  will 
be  one  measure  of  comparison  for  signal  sets  of  different  size,  M,  and  different 
matrices  [ p . j  1 .  The  other  basic  parameter  is  the  ratio  of  bandwidth  occupancy- 
to-data  rate  W/R.  Bandwidth  occupancy  was  shown  above  to  be  n/2T  Hz,  where 
n  is  the  number  of  sub-intervals  of  the  code  signal.  Taking  the  ratio  of  this  to  the 
rate,  it  follows  that 


2  log0  M  '  2K 


In  general,  the  aim  will  be  to  find  classes  of  codes  which  produce  as  low  an 
error  probability  as  possible  for  a  given  K  and  ST^/Nq  ratio,  under  a  constraint 
of  the  W/R  ratio. 


A2.6  Performance  of  Orthogonal  Code  Signals  for  M  —  oo 

Although  orthogonal  signals  are  not  optimal  for  communication  over  the 

Gaussian  channel,  for  large  numbers  of  signals  M  they  nearly  achieve  the  highest 

degree  of  trans -orthogonality  because  the  maximum  correlation  coefficient  can  be 

no  smaller  than  -1/(M  -  1),  which  is  very  nearly  zero  for  large  M.  Also  their 

performance  is  simplest  to  analyze  because  the  correlation  matrix  [p . . ]  =  I  and 

A  2 

consequently  the  expression  for  error  probability  reduces  to 
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(A24) 


Pj2  is  indeptndent  of  which  signal  was  transmitted  since  Vj  is  merely  a  dummy 
variable  of  integration.  The  above  equation  has  been  integrated  numerically  using 
an  IEM  computer  for  M  -  2^  for  K  =  1  through  10  and  for  K  =  15  and  20.  For 
proper  comparison  of  the  results  for  different  values  of  M,  the  independent  v  table 
should  be  ST^/Nq  where  =  T/logg  M.  Figure  A7  presents  results  and  d<  ,on- 
strates  that  for  all  but  very  low  values  of  STg/N  ,  the  error  probability  decreases 
with  the  number  of  signals  in  the  set,  M,  or  the  data-vector  dimension  K. 

It  is  particularly  interesting  to  consider  the  limiting  behavior  of  Pj?  as  M 
approaches  ».  The  similar  results  shown  below,  are  also  true  fcr  noicoherent 
reception  of  M-ary  transmission.  The  optimum  receiver  in  this  case  consists  of 
a  set  of  M-envelope  correlation  detectors  returning  to  the  above  equation 


A2.  Golomb,  S.W.,  editor  (1904)  Digital  Communications  with  Space  Applications, 
Prentice  Hall,  New  York,  KTTY. ,  Chapter  7. 
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Figure  A7.  P„(k)  as  a  Function  of  k,  Showing 
That  for  All  Bat  Very  Low  Values  of  STt>/N„, 

n  _ _ _  urtu.  it.i  ...  i.  U 


Pp  Decreases  With  M  or  k 


M  -  oo  C  '  M 


Li  in  7  ex^-vl/2» 
1  -  00  J  Jilt 


erf  I  +* 


7  exP[-v2/2i  l  im  r  / 

I  M  -  «  J  °  V 1 

-00  L  x 


'2ST^  log2  M 


'2 ST ^  log2  M 


dvj 


Let  STg/N^  (1  and  consider  the  limit  of  the  logarithm  of  the  expression  in  brackets 
of  Eq.  (A25) 
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Taking  M  as  a  continuous  variable  and  using  i  Hospital's  rule  produced  the  above 
limit  as 
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(A27) 


Since  this  is  the  limit  of  the  logarithm  of  the  expression  in  brackets  in  Eq.  (A25), 
the  limit  of  the  expression  itself  is 
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which  implies  of  Eq,  (A25)  that 
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or  in  terms  of  error  probability 
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The  stepwise  limiting  behavior  is  shown  in  Figure  A7.  The  limiting  behavior  is 
particularly  significant  from  an  information  theoretic  vuwpoint.  Shanun's  formula 
for  the  channel  capacity  of  a  channel  of  bandwidth  W  pertuibed  by  additive  Gaussian 
noise  of  uniform  spectral  density  is 
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This  represents  the  maximum  rate  for  error-free  transmission  in  the  limit  as  the 
message  duration  becomes  infinite.  Letting  the  bandwidth  also  approach  infinity 
yields 
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Thus  from  Eqs,  (A30)  and  (A32)  it  follows  that  orthogonal  cedes  achieve  error- 
free  transmission  in  the  limit  as  M  -  »  provided  the  data  rate 
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Of  course,  as  M  -  «>,  the  message  duration  also  approaches  infinity  since 
T  -  Tp  log0  M.  Also  the  bandwidth  approaches  «.  Thus  if  no  bandwidth  constraint 
is  placed  on  the  code  signal  set,  otdiugunal  codes  exhibit  the  best  possible  behavior 
in  the  limit  as  M  approaches  infinity,  since  they  provide  erior-free  transmission 
rates  up  to  channel  capacity. 


A2.7  The  St*<|iieutial  vs.  Hit  Krror  Probability 

The  significant  measure  of  any  communication  system's  performance  depends 
upon  its  use.  If  a  set  of  K-btt  messages  such  as  teletype  or  sampled  data  is  to  be 
sent,  the  sequence-error  probability,  that  is,  the  ptobability  that  a  decoding  error 
was  made  in  at  least  one  bit  in  a  sequence  of  K  bits,  is  the  important  parameter. 

On  the  other  hand,  if  a  sequence  of  independent  bits  is  sent,  the  bit -error  ptobabtl- 
ity  should  be  determined. 

For  orthogonal  coding,  since  ail  emirs  are  equally  probable,  the  expected 
number  of  bits  in  ertoi  when  a  K-bit  i  uded  *otd  has  been  detect e  i  incorrectly  is 


ill 


Thus  the  conditional  probability  that  a  given  bit  is  in  error  when  the  K-bit  word 

K -I 

within  which  it  was  encoded  is  incorrect  if,  2  /(2  -  1). 

Thus  in  terms  of  the  word-error  probability  Pn(K)  for  a  K-bit  orthogonal 

A?  “ 

code  word,  the  bit-error  probability  is 


Pb(K) 


2k-i 


pe(K) 


( 


where  K  =  log9  M.  P.,(k)  is  plotted  in  Figure  A 8  for  various  values  of  k. 


To  conclude,  symbol-error  probability  is  related  to  binary-error  probability 
by  the  formula 


PE(k)  «  (M  -  1)  PE  b 


(A  3  6) 


where  P^  b  is  the  binary-error  probability;  we  now  have 
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or 
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(A  3  8) 


A2.8  M-ary  TmntmUaion  Over  Fading  Channel* 

The  following  situation  is  considered.  A  transmitter  sends  one  of  M  equal 
energy  waveform  over  a  linear,  randomly  time-varying  channel.  These  waveforms 
are  all  frequency  shifts  of  the  same  basic  waveform.  The  channel  is  assumed  to 
have  stationary  statistics  and  stationary  behavior  in  frequency.  The  transmitted 
waveform  will  be  corrupted  in  a  statistically  identical  manner  no  matter  which  of 
the  M  frequency  translates  is  sent.  In  addition  to  the  time  variations  of  the  impulse 
response  of  the  channel  which  cause  a  multiplication  disturbance,  the  signal  is  also 
corrupted  by  additive  white  noise. 

The  receiver  measures  the  received  energy  at  each  of  the  M  frequencies  and 
yields  an  output  symbol  corresponding  to  the  largest  measured  value.  These 
energy  measurements  are  performed  by  correlating  the  received  signal  with  the 
appropriate  frequency  translates  of  a  finite  set  of  orthonormal  functions  and  then 
squaring  and  adding  the  correlator  outputs.  It  will  be  assumed  that  the  frequency 
separation  of  the  M-possible  transmitted  waveforms  is  chosen  sufficiently  large  so 
that  the  signal  energy  appearing  in  one  set  of  correlator  outputs  as  a  result  of 
transmitting  a  pulse  at  a  different  frequency  can  be  ignored. 

Under  these  conditions  Pierce^  showed  that,  for  a  fixed  transmitted  energy 
per  bit,  as  the  alphabet  size  is  allowed  to  become  arbitrarily  large,  the  probabil¬ 
ity  of  correct  reception  depends  only  on  the  noise  density  and  the  probability  dis¬ 
tribution  of  the  energy  pet  bit  contained  in*the  projection  of  the  signal  on  the  sub¬ 
space  spanned  by  the  orthogonal  lunctions  sufficiently  large,  the  error  probability 


A3.  Pierce,  J.N.  (1966)  Ultimat  *  performance  of  M-ary  transmission  on  fading 
channels,  IEEE  Trans,  on  Information  Theory,  Vol.  12,  No.  1. 
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can  be  made  as  small  as  desired  as  long  as  the  expectations  of  total  received 
energy  per  bit  exceeds  (log  2)  times  the  thermal  noise  density.  The  result  is 
somewhat  surprising  in  that  it  is  independent  of  how  badly  the  waveforms  are 
corrupted  by  the  multiplicative  channel  and  of  the  exact  fading  distribution.  It  is 
also  shown  that  this  limiting  capacity  cannot  be  exceeded  by  any  other  type  of  com¬ 
bining  method  if  the  possible  transmitted  symbols  are  equiprobable.  This  is  due 
to  the  fact  that  the  performance  achievable  with  the  pure  energy  addition  is  the 
same  as  <.he  performance  attainable  on  the  infinite-bandwidth  nonfading  channel. 

A3.  CONCLUDING  REMARKS 

There  are  additional  factors  which  have  been  ignored  in  the  present  analysis 
which  tend  to  limit  the  performance  of  actual  communication  systems.  Among 
these  are  distortion  in  the  received  signal  and  inter  symbol  interference  due  to 
nonlinear  delay,  band  limiting,  and  fluctuations  in  the  gain  of  the  medium  that 
couples  the  transmitter  to  the  receiver.  Furthermore,  the  received  signal  is 
processed  in  less  than  ideal  fashion  by  the  detector  due  to  imperfections  in  the 
hardware  and  time  recovery.  It  is  noted  that  all  these  factors  ultimately  manifest 
themselves  at  the  detector  simply  as  a  perturbation  in  the  position  of  the  trans¬ 
mitted  message  point.  As  such,  the  effects  are  similar  to  those  produced  by  the 
noise  and  may  largely  be  compensated  for  by  an  additional  margin  of  signal-to- 
noise  ratio  at  the  detector. 

Conversely,  we  might  say  that  some  fraction  of  the  total  signal-to-noise  ratio 
available  at  the  detector  is  needed  to  compensate  fo1'  effects  of  the  type  listed 
above  which  were  not  accounted  for  in  the  basic  analysis.  Consequently,  only  the 
remaining  fraction  of  the  signal-to-noise  ratio  is  available  for  combating  Gaussian 
noise.  It  is  to  be  expected,  therefore,  that  any  predictions  of  the  probability  of 
error  based  on  estimates  of  the  total  signal-to-noise  ratio  available  at  the  detector 
will  be  unduly  optimistic. 
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Appendix  B 

Preliminary  Ettimatas  of  Error  Rates  in  Asymmetrical,  Doppler- 
Shifted,  Transaurorai,  HF  Propagation  Paths 


This  appendix  discusses  the  influence  of  a  Doppler-shifted  asymmetrical  fad¬ 
ing  spectrum  on  the  binary-error  probabilities  of  Incoherent  FSK  and  differentially 
coherent  phase-reversal  matched-filter  receivers. 

III.  INTRODUCTORY  REMARKS 

Recently,  Gupta®1  has  developed  a  complex  filtering  technique  to  simulate 
Asymmetrical  Doppler-Shifted  Narrow-Band  Gaussian  (ADS-NBG)  channels.  With 
this  technique,  it  became  possible  to  account  for  some  additional  channel  param¬ 
eters  not  considered  before,  such  as  the  non-zero  centroid  and  the  skewness  of 
fading  spectrum  (or  of  the  spectrum  of  the  extended  channel  comprising  the  fading 
channel  and  receiver  processing  filters).  Such  spectral  shapes  are  known  to  occur 
in  troposcatter  channels  (see  references  in  Gupta®*)  and  in  the  HF  transaurorai 
channels,  of  direct  interest  here.  " 


Bl.  Gupta,  A.K.  (1979)  Complex  Envelope  Simulation  of  Scatter  Channels,  pre¬ 
sented  at  the  IEEE  National  Telecommunications  Conference,  Washington, 
DC,  November  27-29. 

B2.  Lomax,  J.B.  (1970)  HF  Propagation  Dispersion,  in  "Phase  and  Frequency 

Instabilities  in  Electromagnetic  Propagation  "  K.  Davies,  Editor,  pp.  497- 
510,  AGARD  Conference  Proceedings,  Number  Thirty-Three,  Techvision 
Services,  Slough,  England. 
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The  reasons  for  this  non-zero  Doppler  shift  and  for  the  asymmetry  in  the 
spectrum  of  the  baseband  process  are: 

1.  Doppler  shift  (see  Figure  Bl).  This  may  be  duo  to  "offset"  in  the  tuning  of 
the  link's  oscillators,  to  time  changes  in  the  refractive/scattering  properties  of 
the  medium,  etc, ; 

2.  Asymmetry  in  the  spectrum  (see  Figure  B2).  This  may  be  due  to  geomet¬ 
rical  asymmetries  in  the  link  (off-great-clrcle-path  superimposition  to  the  main 
path),  asymmetries  in  the  receiver  filters  (RF,  IF  or  baseband),  etc. 

Theoretical  Interpretations  developed  thus  far  account  only  for  the  Doppler 
shift  due  to  the  frequency  offset  in  Ine  local  oscillator  at  the  receiver.  No  analyti¬ 
cal  description  had  been  worked  out  before  in  order  to  account  for  the  shift  due  to 
propagation  medium  and  for  the  asymmetry  in  the  received  baseband  spectrum. 

Figure  B3  shows  a  typical  representation  available  in  the  literature  for  the 

received  spectrum  which  is  clearly  symmetrical.  This  spectrum  is  then  shifted 

by  the  "spectral-shift  unit"  to  simulate  the  frequency  offs  it  in  the  receiver  oscil- 
B 1 

lators.  In  Gupta  an  analytical  representation  of  the  channel  has  been  developed 
that  fully  accounts  for  all  aspects  of  the  two  phenomena  mentioned  above,  thereby 
removing  the  limitations  In  the  present  channel  simulators. 

In  this  section,  we  investigate  the  effects  of  non-zero  centroid  and  skewness 
of  the  fading  spectrum  on  the  binary-error  probabilities  of  incoherent  and  different¬ 
ially  coherent  phase-reversal  matched-filter  receivers  employing  post -detection 
diversity  combining. 

In  this  analysis,  flat  fading  is  assumed.  The  expressions  of  the  binary-error 
probabilities,  derived  here,  are  the  extension  of  the  results  for  the  case  of  the 
Gaussian  fading  correlation  function  of  Bello  and  Nelin.^3  The  general  expres¬ 
sions  for  the  binary-error  probabilities  of  incoherent  and  differentially  coherent 
matched-filter  receivers  ore  derived.  In  Bello  and  Nelin,  the  signal-to-noise 
degradation  (due  to  lading)  curves  are  given  as  a  function  of  fading  bandwidth. 
However,  we  also  obtain  the  signal-to-noise  degradation  expressions  due  to  non¬ 
zero  centroid  and  skewness  of  the  spectrum. 

B4 

Recently  Gupta  has  approximated  an  asymmetrical  Doppler-shifted  fading 
spectrum  by  its  Edgeworth  expansion.  This  expansion  has  the  property  that  the 
second  through  fourth  central  spectral  moments  are  identical  to  that  of  the  spec¬ 
trum  approximated.  If  is  the  centroid  (first  simple  moment)  and  mg,  m^  and 


B3.  Bello,  P.  A.  and  Nelin,  B.  D.  (19G2)  The  influence  of  fading  spectrum  on  the 
binary  error  probabilities  of  incoherent  and  differentially  coherent  matched 
filter  receivers,  IRE  Trans.  Communications  Systems,  Vol.  CS-10,  pp. 
160-168,  June. 

B4,  Gupta,  A.K.  (1979)  Dissertation  abstract,  IEEE  Transaction  on  Information 
Theory,  November,  p.  760.  '  " 
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Figure  Bl.  Doppler-Shifted  Spectrum  in  a  Transauroral  Path  and  an  Example  of 
Computer  Simulation 
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Figure  B2.  Example  of  Skewed  Asymmetrical  Spectrum  in  Transauroral  Path 
and  a  Computer  Simulation 
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Figure  B3.  Typical  Symmetrical  Spectra  Considered  Thus  Far  in 
the  Literature 


m4  are  second  through  fourth  central  spectral  moments,  then  the  Edgeworth  ex¬ 
pansion  of  the  spectrum  S(u)  is  given  by 


(Bl) 


In  Eq.  (Bl)  Hen(y)  is  a  Hermite  polynominal  of  degree  n  defined  by 
Ilen(y)  =  exp  (v2/2)  (-d/dy)n  exp  <-y2/2). 

The  Fourier  inverse  transform  of  S(u)  is  given  by 
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(B2) 
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where,  as  in  probabilistic  analysis,  s^  and  Ex  are  the  skew  coefficient  and  the 
coefficient  of  excess  respectively.  For  simplicity,  we  will  assume  that  Ex  =  0 
and  that  the  last  term  in  the  bracket  of  Eq.  (B2)  is  negligible.  Therefore  Rj(t)  in 
Eq.  (B2)  will  consist  of  the  sum  of  two  terms: 

1.  the  term  involving  the  Doppler-shifted  but  symmetric  spectrum 

2.  the  additive  interference  term  involving  the  asymmetry  of  the  spectrum. 
Therefore  R^('r)  is  approximated  to 

R(t)  =  2a2  exp  ^  jmj  t  -  i  nr^  t^(1  +  aT3]  (B3) 

where  a  =  -j  {s^  m^’  J/6)  denotes  the  magnitude  of  interference. 

Now  we  obtain  the  probablllty-of-orror  expressions  for  the  incoherent  FSK 

and  differentially  coherent  phase-reversal  matched-filter  receivers.  It  is  to  be 

B3 

noted  that  in  Bello  and  Nelin,  m  j  and  sk  are  assumed  to  be  zero. 


1)2.  ERROR  PROBABILITIES  FOR  INCOHERENT  MATCHED  FILTER  RECEIVERS 


A  simplified  biock  diagram  of  an  incoherent  matched-filter  receiver  is  shown 

In  Figure  B4.  Two  matched  filters  are  Indicated,  each  matched  to  one  cf  the  two 

possible  waveforms  that  are  successively  repeated  to  form  the  transmitted  binary 

communication  signal.  The  outputs  of  the  matched  filters  are  envelope-detected 

and  subtracted.  Diversity-combining  is  performed  by  summing  the  outputs  for  all 

the  diversity  receivers.  It  is  assumed  that  bit  synchronism  exists.  Based  upon 

the  stationary  complex-valued  Gaussian  fading  statistics,  the  error  probability  for 

a  received  binary  waveform  occupying  the  interval  of  bit  duration  are  calculated  in 
B3 

Bello  and  Nelin.  For  the  case  of  FSK  transmission,  the  transmitted  waveforms 
are  given  by 


(B4) 


where  the  integer  n  is  equal  to  the  frequency  separation  between  the  mark  and  space 
frequencies  normalized  with  respect  to  the  data  rate  1/T. 


Input  Signal 


Matched 
filter  1 

recei vers 

hH 

Envelope 

detector 

1 _  1 

_ _ a 

7 

Matched 
filter  2 


Envelope 

detector 


From  other 


Subt 

ract  ■ 

5 

Sample 

From  other 
recei  vers 


Figure  B4.  Block  Diagram  of  an  Incoherent  Matched-Filter  Receiver 


For  the  waveforms  given  in  Eq.  (B4),  the  required  m  and  Xjig  functions, 
are  defined  in  Bello  and  NelinE3  as 
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T 

Ers  =  j  Ss(t)  Sr(t)  dt  ,  where  Nq  is  the  rpecira)  density  of  additive  noise, 
o 

and  S(t)  could  be  either  Sj(t)  or  SQ(t),  are  given  by 
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9  [  -2jnir 

•  °-<T 

2  [  ^SiSIr' 

-Sr  L‘-°  J  •  -T<T<0  • 


In  Eq.  (B6),  the  superscript  "o"  denotes  S(t)  =  SQ(t),  o  <  t  <  T,  "1"  denotes 
S(t)  =  Sj(t),  o  <  t  <  T  and 


E  -  i  /  |So(t)|2  dt  =  1  f  |sl(t)l2  dt  .  0 

o  o 

2 

Substituting  R(t)  in  Eq.  (B3)  in  m  Eq.  (B5),  one  obtains  (letting p  =  cf  E/Nq) 

-(l/2)m2T2 

9  9  ,  i  2  cos  (m.T)  e  -2 

m°  (t)  =  16  E“  a  ~  +  - =r“T - 

oo  2p  21  nio 1 


\  m2i  /  0  \  m2*  /  0  rrij  replaced  by 

,  8  E2  o2  -l°jmlm2  f  2jmi  .  T  T 

+  — ™ -  a - 2 - +  *i  "  M 

1  mgT  nij  replaced  by  -mj 


+  l2  ”  *2 


m,  replaced  by  -m, 


In  Eq.  (B8) 
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l  ,  2nir\ 
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where 


T3  - 


2m 2  -  m^ 
m?, 


2m2  -  m\  +  T2m2  +  jm^T 


mr 


* 


•  exp 


T 


Similarly  one  can  obtain 


m 


1 

oo 


m 


11 


n  replaced  by 
-  n 
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mj.  =  m° 
11  oo 


(B13) 
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m,  replaced  by 

•  (mi  -~t) 


113.  PROBABILITY  OF  ERROR  FORMULATION  FOR  INCOHERENT  RECEPTION 

B3 

To  compute  the  probability  of  error,  one  requires  parameter  r  as 


(B 14 ) 


2(>nil  -  m<K>) 


*/ [m, ,  -  m  )2  +  4  (m.  ,m  -  | m ,  I  -  (m. .  -  m  ) 
v  11  oo  11  ui  lo  11  oo 


(131 5) 
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Depending  upon  the  superscripts  on  mr  s(r,  s  =  o,  1),  one  obtains  r  ,  r^.  The 
significance  of  rQ  and  r^  will  be  outlined  below. 

The  probability  of  error  expression  for  zero-fading  bandwidth  and  square-law 
combined  diversity  was  derived  by  Pierce.  if  is  the  probability  that  a  "o"  is 
printed  given  that  a  11 1"  is  sent  and  pQ  is  the  probability  that  a  "  1"  is  printed  given 
that  a  "o"  is  s«nt,  then  for  zero-fading  bandwidth  pQ  =  p1  =  p  where 


M-l 


m=o 


(2M  -  l)l(-l)m _  ,,  ,  rM-m 

(M  -  HI  (M  -  1  -  mTs  ml  (M  +  m)  p) 


(BIG) 


In  I2q.  (BIG),  M  is  tne  order  of  diversity  and  p  is  the  ratio  of  the  slgnai  power  in 
data  channel  and  the  noise  power  in  data  channel.  The  noise  power  in  the  data 
channel  is  defined  as  the  noise  power  in  a  rectangular  bandpass  filter  of  unity  gain 
and  bandwidth  l/T.  At  a  large  slgnal-to-noise  ratio,  SNR,  the  leading  terms  pre¬ 
dominate  and  we  have,  to  a  good  approximation 


P 


(2M  -  1)! 
(M  -  l)!  Ml 


p  •»  1 


(B17) 


B3 

For  ».hc  case  of  general  fading  and  general  binary  signals,  one  may  deter¬ 
mine  Pj  by  using  an  "equivalent"  slgnal-to-noise  ratio  Pj  in  the  zero-fading  band¬ 
width  expression  for  pj(Eqs.  (BIG)  or  (B17)),  This  equivalent  SNR  is  given  by 


p.  ^  r 


(BIG) 


Similarly  pQ  for  general  fading  can  be  determined  by  using  an  equivalent  SNR  pQ,  in 
the  slow  fading  expression  for  Pj  where 


po 


1 


r 

o 


t-  r 

o 


(B19) 


In  general  pQ  *  Pj.  Thus  in  the  presence  of  an  asymmetrical  spectrum  and/or 
a  non-zero  centroid,  asymmetric  binary  channel  may  result. 


B5.  Pierce,  J.N.  (1958)  Theoretical  Diversity  Improvement  in  Frequency-Shift 
Keying,  Proc.  IRE,  Vol.  45,  pp.  903-910,  May. 
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B4,  ERROR  PROBABILITIES  FOR  A  PHASE  REVERSAL  KINEPLEX  SYSTEM 


A  simplified  block  diagram  of  the  differentially  coherent  matched-filter 
receiver  is  shown  in  Figure  B5.  It  is  assumed  that  only  two  possible  pulses  S ^ (t>, 

«  SQ(t)  are  transmitted  as  with  the  incoherent  matched-filter  receiver.  However, 

the  transmitted  bit  is  encoded  into  the  change  or  lack  of  change  of  successive  adja¬ 
cent  pulses.  Thus  a  transmission  of  the  pairs  Sj(t)  or  SQ(t>,  Sj(t)  denotes  a  "o", 

*  say,  while  the  pairs  Sj(t),  S^(t)  and  SQ(t),  SQ(t)  denote  a  "  1" .  The  receiver  con¬ 

sists  in  part  of  two  separate  coherent  matched-filter  receivers  with  filters  matched 
to  waveforms  Sj  and  SQ  but  differing  by  an  input  delay  of  one  band  duration.  Thus, 
the  two  receivers  are  operating  upon  adjacent  bands,  The  actual  detector  output  is 
obtained  by  using  the  output  of  one  coherent  receiver  as  a  reference  for  the  other. 

In  the  special  case  when  he  pulses  S^t),  S^t)  arc  ±  1  over  the  band  duration,  the 
differentially  coherent  system  specializes  to  the  differentially  coherent  phase- 
reversal  Kineplex  system. 


Figure  B5.  Block  Diagram  of  a  Differentially  Coherent  Matched-Filter  Receiver 


10' 


In  the  following  discussion,  we  will  compute  the  expression  of  the  probability 
of  error  for  DPSK  or  phase-reversal  Kineplex  system,  Due  to  antipodal  nature  of 
binary  waveforms 


S^t) 


and  SQ(t)  =  -Sj(t)  for  o  <  t  <  T  , 


(B20) 


binary  symmetric  operation  will  exist  {that  is,  pQ  =  p^)  for  all  fading  correlation 

functions.  As  in  the  FSK  Incoherent  case,  the  following  results  are  the  extension 

of  the  results  given  in  Bello  and  Nolin.  The  reader  is  referred  to  Bello  and 
B3 

Nelin  for  a  detailed  analysis. 

For  the  case  of  general  fading  and  general  binary  signals,  pQ  and  pt  can  be 
determined,  similar  to  the  FSK  case,  by  using  "equivalent"  SNR  in  Eqs.  (BIO)  or 
(B17).  The  error-probability  expression  in  12q.  (BIG)  was  derived  for  noncoherent 
FSK  reception  with  "square  law"  diversity  combining  and  slow  fading  (zero-fading 
bandwidth).  The  equivalent  SNH  is  glven*^  ns 


'll 


2(mlo  +  mio) 


V^o^lo*2*  4(mllmoo  '  !,nJ2)  '  (mlo  +  mlo) 


(B21) 


where 


moosmll  s32e2°2 


l  +  L 

P  I 


-(l/2)m.,T“ 

2  cos  (mjT)e  -2  /  jnij 

m0T  +  I  1  - 


*f 


nij  replaced 
by  -  nij 


■f 


a 

T 


-10j|in^m9  +  2jnij 

mjT 


+  I, 


+  h  -  h 

nij  replaced 
by  -  nij 


ntj  replaced 
by  -  nij 


(B22) 
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and  letting  =  rnj  +  jm2  T 


,  .  32E 2  a2 

*lo - T -  e!£P 


<?  (jmjT  -  I, 


2  cos  (m'jT)  e 


-(l/2)m0T2 
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T4  +  T3  + 


6m' 
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m. 
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6T 
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mr 
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m. 
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nij  replaced  \ 
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(B23) 


To  recapitulate,  the  probability  of  error  for  a  skewed  and  Doppler-shifted 
channel  (incoherent  FSK  modulation),  can  be  computed  as  follows: 

1st  Step— Compute  the  Doppler  shift  of  the  received  spectrum.  This  is  called 
m.  (rad  sec”1): 

«  rul 

2  Step— Compute  the  parameter  a,  as  defined  by 


a  = 


-  3 


V"2*5 


where 

m,,  is  the  fading  rate  (-  Doppler  spread  of  the  received  spectrum,  in  rad/sec) 
s^  is  the  skewness  of  the  spectrum,  as  defined  in  Eq.  (22),  Section  2,4,2. 

3 1(1  Stop -Compute  the  quantities  m  (where  r  and  s  acquire  the  values  0  and 

1)  by  using  Eqs.  (BIO)  through  (B14).  (N’ote  that  m®  g  refers  to  bit  "0"  and  m* 
refers  to  bit  "  l"). 

4!b  Step -Compute  the  function  r  from  Eq.  (B15). 

5* Step -Compute  p5  and  pQ  from  Eqs.  (B18)  and  (B19)  (p  refers  to  bit  "0"  and 
p  j  refers  to  bit  "  1"). 

Qlb  Step— Compute  p}  and  p  (equivalent  binary  probability  of  error)  from 
Eqs.  (BIG)  or  (B17)  (again,  Pq  refers  to  bit  "0"  and  Pj  refers  to  bit  "1"). 

Similar  to  what  has  been  just  described  for  the  incoherent  FSK  modulation 
case,  the  probability  of  error  can  be  computed  for  the  case  of  a  coherent-differen¬ 
tial  Kineplex  system.  Since,  in  this  case,  pQ  and  Pj  for  all  shapes  of  the  fading 

spectrum,  the  computation  of  either  one  of  them  is  sufficient.  First,  m,,  ,n, 

it  lo 

and  moo  are  computed  by  using  the  Eqs.  (B22)  and  (B23)  as  a  function  of  the 

spectral  parameters.  Then,  p  ^ j  is  computed  from  Eq.  (B21).  Since  p^  can  be  * 

thought  of  as  an  "equivalent  signal-to-noise  ratio,"  the  probability  of  error  can  be 
computed  from  the  Pierce  formula,  given  in  Eq.  (BIG), 
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